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ABSTRACT 

An attempt is made at developing a radar parameter sufficiently 
correlated with hailfall at the surface that it may be used for cloud-seeding 
evaluation purposes. 

The parameter is volume-averaged radar reflectivity 7, applied 
to low-level PP! scans. The parameter 7 is tested on 9 Alberta hailstorm 
case studies from 1975 and 1976, using Alberta Hail Project 10.4 cm wavelength 
radar data. Cumulative distribution curves for n according to maximum hailsize 
reported at the surface are obtained. The problems associated with 7 are com- 
parable with those of equivalent radar reflectivity for hailstorms, although, 1 
appears to be sensitive to, and a good indicator of, overall storm intensity. The 
quality of the 1975 and 1976 radar data is discussed. 

A spectral analysis of the 7 time series for the storms reveals possibil- 
ities of preferred cycles. Fluctuations in n imply a storm concept consisting of 
a large circulation cell with periodicity of approximately 30 minutes, upon which 
is superimposed small-scale perturbations with possible periodicities of approximately 
10 and 15 minutes. 

The analysis of changes in 7 with respect to cloud-seeding events, 
suggests a positive cloud~seeding effect in some storms, however, unequivocal 


seeding effects were not detected. 
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CHAPTER 1 
INTRODUCTION 
1.1 . Background on Hail Suppression 

The hailstorm poses a serious problem for agriculture in the Province 
of Alberta with annual crop losses averaging $40 to $50 million (Renick, 1975). 
Agricultural claims for hail damage are received on an average of 50 days each 
season (1 June to 10 September) in Alberta (Summers and Wojtiw, 1971), with 
50 percent of the years' damage occurring on 4 days, and 80 percent of the 
damage on the worst 12 days (Wojtiw, 1975a). 

There are hypotheses which suggest that hail may be suppressed by 
injecting (seeding) ice nucleating agents such as silver iodide (Ag!) into growing 
convective clouds. The claim is made that the growth of large damaging hail 
is inhibited through the introduction of high concentrations of ice nuclei, all 
competing for a finite supply of supercooled water. Considering the agricultural 
losses, even partially successful hail suppression programs would be economically 
beneficial; thus numerous hail suppression experiments have been carried out in 
North America and abroad. Claims of decreases in hail damage have been made 
in a few cases, but due to the lack of significant evidence in support of such claims, 
no general agreement exists among the scientific community about their validity 
(Hitschfeld, 1974; Neyman, 1977). 

The Alberta Hail Project (AHP), under the auspices of the Alberta 
Weather Modification Board (AWMB) with participation by the Alberta Research 
Council (ARC), is among the most significant hail suppression programs currently 


in progress. The AWMB established under Alberta Agriculture in 1974, was formed 
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in response to the leconaaeitlinns of the Special Legislative Committee on Crop 
Insurance and Weather Modification to the Legislative Assembly of Alberta 
(March 1973). It was granted a five year mandate to administer hail suppression 
and hail research programs with continuing evaluation to determine the scientific 
and economic feasibility of hail suppression using current technology. An area 
of 47360 km* (18,500 mi*) centered on the radar site at the Red Deer Industrial 
Airport was designated for the project (Figure 1.1). Although the AHP cloud 
seeding experiment began in 1974, its organization and many of its researches 
are extensions of the former Alberta Hail Studies (ALHAS) project which originated 
in 1956. The ALHAS project was devoted entirely to basic hailstorm research 
and began experimental seeding of storms with Agl for hail damage reduction in 
1970. 
1.2 1975 - 76 Alberta Hail Project Seeding Techniques 
The AHP cloud seeding program runs from the middle of June to the 
middle of September each year. The cloud seeding is conducted by seven turbo- 
charged twin-engined aircraft either from cloud top, utilizing dropable pyrotechnic 
flares, or from cloud base with pyrotechnic flares mounted on the aircraft wings. 
The cloud~top seeding method was first developed by Simpson et al. 
(1970) in Florida, and adapted to Alberta hailstorms by Summers et al. (1972). 
It is closely related to the Alberta multi-cell storm concept as described by Douglas 
and Hitschfeld (1959), Chisholm and Renick (1972), and Marwitz (1972). 
In 1974 and 1975, seeding was confined to growing cloud towers and 


the main storm-inflow area at cloud base (Figure 1.2). Seeding rates varied but, 
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Fig. 1.1: Map showing the area of the Alberta Hail Project. 
(From Deibert (ed.), 1977) 
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Fig. 1.2: Cloud seeding techniques for 1975. (From Deibert (ed.), 1976). 
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Fig. 1.3: Cloud seeding techniques for 1976. (From Barge et al., 1976). 
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in general, cloud-top flares were dropped from approximately the -1 2°C level 

at 5 second intervals (approximately 300 m). Base flares were ignited once every 
7 minutes unless the updraft at cloud base exceeded 5 ms ~, when 5 base flares 
were burned simultaneously. Both methods were used when possible with a higher 
priority given to cloud-top seeding; however, storm type and amount of cloud 
cover often determined which seeding technique was used. 

Seeding procedures were changed slightly in 1976 (Figure 1.3), as 
a result of radar studies performed at Alberta Research Council (Barge et al., 
1976). It is now felt that seeding the smaller growing cloud towers at the edge 
of the storm more heavily at an earlier stage will give better results. Seeding 
was performed by aircraft flying a "figure 8" flight pattern above the new growth 
region at ae Rea the -10°C level and dropping Agl pyrotechnic flares 
approximately 2 km apart. If the shelf cloud was obscured by mid-level cloud, 
seeding was confined to the weaker updrafts around the edge of the storm at cloud 
base. 

In 1975 the cloud-top seeding flares were 50 gm Agl, WM-105 
formulation from Olin corporation. These flares fall approximately 1800 m during 
their one-minute burn time. The cloud-base seeding flares were 70 gm Agl, 
SR-14 fusees (Nuclei Engineering Incorporated) which burn for 7 minutes. 

In 1976 the cloud-seeding flare types included those from 1975 plus 
70 gm Agl, TB] cloud-top flares (Nuclei Engineering Inc.) which have a burn time 
of 35 seconds after an 11 second delay, and TBI cloud base fusees (Nuclei 


Engineering Inc.) which contain 150 gm of Agl with a burntime of 4 minutes. 
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Cloud-seeding operations both years were carried out according to the 
following criteria (Deibert ed., 1976). The procedure was to declare a day 
"experimental" at the first appearance of a 35 dBZ radar echo inside the Tae 
area. Then the next member in a random series of envelopes was opened and the 
"seed" or "no-seed" instruction contained therein was carried out. If the in- 
struction was "seed", all storms with radar reflectivities greater than 35 dBZ oc- 
curring within the project area were seeded for the remainder of the day. If the 
instruction was "no-seed", storms in the northern area were not seeded unless they 
were moving southward and expected to enter the southern area, in which case 
the no-seed decision was overruled and the storms seeded. 

A buffer zone was established between the northern and southern areas 
to facilitate earlier seeding of southward-moving storms close to the southern area 
border. Storms that formed in this zone and were forecast or observed to move 
south or east were seeded as though they were in the southern area. If the storm 
was forecast to track northward, the zone was considered as part of the northern 
area and randomization procedures were used. 

On 23 July, 1976 the AWMB authorized the project management to 
extend the limits of the northern project area to 176 km (110 mi) towards the north- 


' grea, and enable earlier 


west in order to include what is called a storm "breeding' 
seeding of storms moving into the area from that direction. 
In the southern area, the 35 dBZ radar reflectivity criterion for initiat- 


ing seeding of all storms was used. However, when the forecast called for large 


hail, seeding often commenced prior to the first appearance of a 35 dBZ echo. 
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1.3 Alberta Hail Project Evaluation 

Evaluation of the AHP to date has been based primarily on statistical 
tests on hailfall measurements from "seed" and "no-seed" hailstorms from the 
randomized experiment (Wong, 1975, Wong, 1976). However, due to the high 
degree of variability in space and time of the output of hail from a storm, large 
numbers of seeded and unseeded cases must be observed in order to detect any 
differences between the two with adequate statistical confidence. It appears that 
unless the seeding is very effective, conventional statistical analyses require 
periods of time which may be unacceptable in many places (Schickendanz and 
Changnon, 1970; Goyer, 1975). The probability of evaluating statistically the 
efficacy of cloud-seeding by the AHP using present techniques was reduced after 
the decision was made e the AWMB to seed all potential hailstorms in both the 
northern and southern sectors in 1977, thus ending the randomized portion of the 
experiment after only 3 years of operation. 

The field of hail suppression is in need of some fresh evaluation methods 
which are able to abridge the excessive time requirement imposed by a conventional 
statistical evaluation of hailfall measurements. 

1.4 Radar and Evaluation 

Weather radar has long been considered essential in weather-modifica- 
tion experiments because of its value as a short-term forecasting tool and as a basis 
for directing the field operations; but only recently has it been considered as a 
source of data to evaluate the seeding effects. The major uses of calibrated radar 


data, either in real time, or in subsequent analysis, as given by Changnon and 
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Morgan (1976), include: 


LE 


Ze 


Delineating developing hailstorms for directing seeding activities. 
Counting hailstorms in seeded and nonseeded areas as part of an 
evaluation. 

Monitoring various physical changes in storms before and after 
seeding for evaluation. 

Determining the surface area covered by hail for evaluation and 
directing post-storm surveys. 

Ascertaining the echo volumes containing hail for direction of seed- 
ing activities. 


The usefulness of radar as an evaluation tool is that it is able to scan 


a large volume in space in a relatively short period of time, and is able to provide 


quantitative information concerning echo structures and meteorological processes. 


An objective evaluation must be based on a method which is able to monitor the 


hail-producing capacity of a hail storm prior to, during, and after cloud seeding 


not just the hailstorm results such as total hail mass, or the effects such as total 


crop damage, which are merely functions of this hail-producing capacity (Goyer, 


private communication). Radar may be able to accomplish this. 
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GHAPTER 2 
THEORY 
2.1 Introduction 

The central problem described here is to find a radar parameter, suf- 
ficiently correlated with hailfall at the surface, to warrant its use as an evalua~ 
tion parameter. 

It is particularly difficult to find clear relations between radar para- 
meters and precipitation type and rate for thunderstorms since they contain un- 
known and variable combinations of rain with wet, dry, or slushy hail of various 
sizes and concentrations. In spite of these difficulties, useful relations between 
radar parameters and hail have been noted by Geotis (1963), Barge (1972), 
Changnon and Morgan (1976), and Eccles (1976). 
2eceeine erecr nical Radar Equation 

The meteorological radar equation relates the precipitation particles 
which comprise the illuminated "target" to the echo power received. Probert- 
Jones (1962), showed that the meteorological radar equation is obtained by sub- 


stituting the "radar reflectivity" , of small spherical drops, as given by (2.2.1): 
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where 6, are the radar cross-sections of all the individual scatterers, ™ = 3.14159, 
€ is the dielectric constant, r is the particle radius and 4 is the wavelength of the 
radar being used; into an appropriate form of the radar equation for a distributed 


target given by (2.2.2): 
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where c is the speed of light (3 x 10° msea) ee 


, is the average returned power, P; 


is the peak transmitted power, T is the pulse duration, G is the antenna gain, © 
and % are the azimuth and elevation beamwidths respectively, and R is the distance 


from the radar to the target. The result is the meteorological radar equation (2.2.3): 
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The equation is often written in terms of the radar reflectivity factor Z. By defi- 


nifion: 


MS 2 
Np OS ee ee (2.2.4) 
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where D; are the diameters of all the individual scatterers comprising the target, 
and N; the number concentration. The meteorological radar equation in terms 


of the radar reflectivity factor is: 
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According to Smith (1970), the most important assumptions underlying 


the meteorological radar equation are: 
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The scatterers are homogeneous dielectric spheres, with diameters 
small compared to the wavelength, allowing Rayleigh scattering 
theory to be used. This applies to water spheres with D; <0.07 A 
(Gunn and East, 1954), and to ice spheres with D; <0.16 > (Ryde, 
1946). 

The illuminated volume is filled with scatterers. 

The radar reflectivity factor Z is uniform throughout the illuminated 
volume. 

The scatterers all have the same dielectric constant, that is, they 
are all water drops or all ice spheres. 

Multiple scattering is negligible. 

The main lobe of the antenna beam pattern can be described adequately 
by the Gaussian function. 

Microwave attenuation along the propagation path between the radar 
and the target is negligible. 

The incident and back-scattered waves are linearly polarized. 


Where one or more of the assumptions break down (2.2.4) can be used 


to define an "equivalent radar reflectivity factor" Z,. Z, is sometimes referred 


e 


to as a measure of "what the radar sees", and can be considered to be a parameter 


which is calculated from the radar measurements. For hailstorms (2.2.5) should 


more appropriately be written as: 
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2.3 The Evaluation Parameter 

According to the hypothesis being tested by the AHP, the size distri- 
bution of hailstones will be shifted towards smaller sizes and their number con- 
centration will be increased as a result of effective seeding. This stems from the 
hypothesis that the large, damaging hailstones occur because of a natural defi- 
ciency of ice nuclei, and that the injection of suitable numbers of ice nuclei 
(Agl) into the storm promotes the formation of a larger number of smaller, less 
damaging hailstones, many of which could melt during their fall to the surface. 

If the hypothesis is correct, what changes in the radar reflectivity 
are expected? From its definition in (2.2.4), the radar reflectivity factor varies 
as the sixth power of the target diameter and only as the first power of the number 
concentration of the eontiarers: Eccles (1976), found that for hail below the 


0°C level the relation is: 


eee) Deo 


and above the 0°C level: 


These are very close to the Rayleigh relation. 
As an example, assume that a volume in space contains a hailstone 
concentration of .01 m7*, with each hailstone having a diameter of 20 mm. The 


radar reflectivity factor is: 
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= 6.40 x 10° mm° m7~ 


58.06 dBZ 


According to English (1973) the natural concentration would have to be increased 
by at least a factor of 100 in order for depletion to become significant. Assuming 
that the supply of supercooled water is fixed, the new diameter of the hailstones 


will be D, where: 


The new radar reflectivity factor will be: 


(1.0 m7*) (4.31 mm)° 
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The radar reflectivity factor has been reduced by a factor of 100 (20 dBZ). There- 
fore, effective seeding should lead to a decrease in the reflectivity of the volume 
which contains the hailstones within a storm. 

Nonetheless, reflectivity alone has been shown to be a rather in- 
consistent indicator of hail (Barge, 1972; Changnon, 1972; Martner and Dye, 
1978). However, it has been found by Chisholm (1968) and Changnon and Morgan 
(1976), that there is a good correlation between high surface reflectivity areas 
on a given day and the area of hail at the surface, and that hail shafts near the 


ground generally exist in steep reflectivity gradients. 
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With this in mind, the parameter chosen for study was radar reflect- 
ivity "1, averaged over the volume of the storm echo near the surface. The mod- 
ifier "near the surface" in this case means echoes returned fromthe 1° to 2° ele- 
vation angles of the radar antenna, but not including ground clutter. Table 2.1 
gives the height above ground of the radar beam for a given range at low eleva- 
tion angles. The reflectivity values received from the low levels can be assumed 
to come from precipitation scatterers which actually reach the ground. In this 
way the parameter utilizes the good correlation between high surface reflectivity 
areas and hailfall areas at the ground. Volume averaging is used in order to 
take into account the fact that the radar beam diverges with increasing range and 
that the reflectivity values are determined over larger volumes as the target - 
antenna distance ince: 

A useful expression involving 7 can be found in (2.2.3). This equation 
expresses the average received power P; in terms of the average radar cross section 
or reflectivity 1 of the array of scatterers distributed throughout the volume illumi- 
nated by the radar beam. This form is preferred because all terms are either known or 
easily measured, and all uncertainties are associated with the radar cross section only. 

According to Smith (1970), the radar cross section depends on the 
following characteristics: 

1. Shape 

2. Size 

3. Dielectric constant 


4. Viewing aspect 
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Range 


100 


Table 2.1: Radar Beam Height (meters above ground) 


km 


o2 


48 


64 


80 


96 


HZ 


128 


144 


160 


From equation h = 


Elevation Angle (degrees) 


le io ze 
619 ize {| eke 
973 1811 2647 
1358 2474 | 3590 
1773 3168 | 4562 
2218 3892 | 5565 
2693 4646. | 6597 
3198 5430 | 7660 
3733 6244 | 8752 
4298 7088 | 9875 


R(Riee2ireesimac) 
r, tvr? +R(R +24, sin 8) 


elevation in degrees 


range in kilometers 
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f earth's radius (mean radius of earth is 6371 km). 


Beam Width 
(meters) 

642 

964 
1285 
1606 
O27, 
2248 
2569 
289 | 
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lt should be noted that the radar cross section bears no universal relationship to 
the actual geometrical area of the target even though its dimensions are (mm? m-*). 
The proposed evaluation parameter therefore is based on ™ the radar 


reflectivity, where: 
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Applied to a storm situation the evaluation parameter is volume-averaged radar 


reflectivity 7 , defined as: 


v 


Jnay 


f dV 


where volume averaging is carried out over the region in space covered by the 

:e) ° * . . * . . . . 

1” to 2: elevation scans. It is proposed that 7 is sensitive to reflectivity changes 
and is able to indicate the presence of steep reflectivity gradients, which as men~ 
tioned, have also been associated with hail. This method presents a frequent 


—_ 


measurement of the hail producing capacity of the storm. Changes in 7 with time 


are analyzed in an effort to detect cloud seeding effects. 
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CHAPTER 3 
THE RADAR FACILITIES 
3.1 The Alberta Hail Project Radar System 
Figure 3.1 depicts the general aspects of the AHP radar system. All 
of the radar data for this thesis were obtained with the AHP S-band (10.4 cm 
wavelength) radar. The characteristics of the S-band radar are as follows (Barge, 


1974; Leung, 1977): 


frequency 2.880 GHZ 

pulse repetition frequency 480 s~* 

pulse duration 1.75 us 

peak power (nominal) 200 kw 

beam dimensions 1.15° (all planes) 

antenna mes 43.2 dB (at pedestal) 

antenna rotation rate 8 rpm 

elevation program Spiral scan: ie per revolution 


to 8° in 1.5 min, or 20° in 

3.0 min. 
The S-band radar is equipped with a variable polarization antenna and a dual 
channel receiver, and is capable of transmitting elliptically polarized radia- 
tion at any chosen axial ratio and orientation. The received radiation is resolved 
by the antenna microwave circuits into its main and orthogonal components which 
then enter separate receiving channels. Only the main component was used in 
this study. 


In addition, a C-band (5.4 cm wavelength) radar is used in a surveillance 
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mode, or when the S-band is turned off for servicing or because of strong winds. 
Finally an X-band (3.2 cm wavelength) radar is used for tracking aircraft. 
3.2 The Alberta Hail Project Radar - Computer System 

A PDP 11/50 computer system is used to record radar data on magnetic 
tape. Using calibration information, the numbers stored on tape are converted 
into received power (Barge, Humphries, and Johnson; 1976). The parameter tested 
here was applied to computer-recorded radar data, therefore it is important to 
know the characteristics of the data. 

For each radar pulse the return signal is detected and logarithmically 
amplified. It is then averaged in 7 Us intervals or "bins" which correspond to a 
length of 1.05 km. The height and width of the bins are determined by the radar 
beam dimensions at a given range. A bin volume thus represents a volume in space 
for which a certain radar reflectivity has been measured. Figure 3.2 depicts the 
bin volume and Figure 3.3 shows a storm echo as it is recorded on magnetic tape. 


The equation for calculating 1 from the stored data is: 


a D 2 yy, 
—_ bins : Ae ied 
iq = 
x Ve 
bins 
On oe 
rae oe eae 


and the bin volume is given as 


V; = R 68 AR AR = 1.05 km 
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The summation is over the total number of bins which comprise the echo. 
3.3 Error Analysis 

The AHP radar is considered to be a stable system and its specifica- 
tions fluctuate very little (Humphries, Barge and Johnson, 1977). Since the 
analysis deals mostly with changes in 7, the primary concern is with "relative 
errors". Although the value of 7 calculated for the AHP radar system (or most 
meteorological radars, for that matter) may have an absolute error of 1 3 dB 
(50%), the relative error is considerably less. The absolute error would ap- 
propriately be adopted when comparing 7 values with those obtained by a radar 
system elsewhere. 

The main source of relative error is "signal fluctuation "due to a 
fluctuating echo with she errors assumed to be constant. From theoretical cor- 
siderations Marshall (1971) showed that if, on a plan-position indicator (PPI) 
display, every independent datum was displayed separately with signal level re- 


ported accurately, the resolution in position would be as good as theoretically 


possible and the target intensity level would be known everywhere with a standard 


deviation 5 = 5.57 dB. If the data are combined and averaged in groups of j 
independent data, the standard deviation is related to the number of independent 


data by: 


D0) ab 
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The AHP S-band radar has 4 independent samples per range bin and 
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10 pulses per degree azimuth (Barge, Humphries, and Johnson, 1976). The com- 
puter averages 10 pulses for each range bin; consequently 40 samples of the re- 

turned power are averaged. Not all of the samples are independent, however, 

since decorrelation in azimuth occurs over one beam-width. Moreover pulse-to- 
pulse frequency shifts of the magnetron and "shuffling" of precipitation also lead 
to further decorrelation (Atlas, 1964). The number of independent samples con- 
tributing to each average is assumed to be 16 and therefore, the target intensity 


level would have a standard deviation of: 


oo ES ee 


16 


Since this analysis involves averaging over a number of bins, for a 
given echo the appropriate standard deviation for 1 may in Sti be decreased 
by the factor 1//m where m is the total number of bins. Figure 3.4 shows the re- 
lation between the relative error involved in 1 and the number of bins in the 
echo. Since only well-developed storms were seeded, the number of bins usually 
exceeded 100. Therefore during the times of interest a relative error of less than 
* 0.139 dB (3%) may be assumed. Adding the relative error due to transmitter 
fluctuations, which may be assumed to be no more than * .5 dB (12%), the relative 


error involved with 1 is not more than * 0.639 dB (15%). 
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Relative error (percent) in n calculation as a function of the 
number of bins in a given PPI echo. 
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CHAPTER 4 
DATA ANALYSIS 
4.1 Analysis of Digital Radar Data 
The computer was programmed to perform the calculation of 7 (see 

Appendix C for program). However, before numerical operations were carried 
out the digitized data had to be properly prepared. This operation was accom- 
plished with the AHP interactive computer display facility. Briefly, the inter- 
active display facility provides the following capabilities (Ramsden et al., 1976): 

1. The data are displayed in the conventional PPI format. 

2. The data analyst is provided with a means to define a portion of the 
PPI, both for display purposes and for extraction of data for sub- 
sequent numerical processing. 

3. Any portion of the data may be displayed at any chosen scale in any 
desired degree of detail up to the maximum resolution provided by 
the data. 

4, Permanent copies of any display are available immediately. 

The interactive display facility is operated from a cathode-ray-tube 
(CRT) graphical display terminal, and uses a CRT screen both for input of instruc- 
tions and for output of graphical displays. Figure 4.1 shows the basic display which 
consists of the PPI display, the antenna elevation display, and the command list. 
The command list is displayed at all times and input of a command is accomplished 
simply by positioning a cross-hair cursor over a desired command in the list and 


pressing any key on the keyboard. 


Initially the display consists of a complete PPI. However, the analyst 
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DISPLAY OPTIONS: 
FIDUCIALS OFF ON 
WINDOWS OFF ON 
LANDMARKS OFF ON 
CONTOURS OFF ON 
GRND REF SET CLR 
OUTPUT OPTIONS: 
HARD COPY OFF ON 


WINDOW OFF ON 
DATA OFF ON 
COMMENT 


DATA TYPE DBM DBZ 


WINDOW OUT IN 12345 = 
VIEWPORTS 1" 4 9 


GO TO TIME 
BEREAt GRREU BNEXT 
SKIP BCK TO 
SKIP FWD TO 
UN 
8. S. 10. 1s. 20. 


BEGR ELEVATION 
19°2e6'36 1 3 DEGR 


Fig. 4.1: The basic radar display using the interactive computer display 
facility. (From Ramsden et al., 1976). 
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Fig. 4.2: The display shown in Fig. 4.1 after a portion of the PPI has 
been outlined for further detailed study. Note the cross-hair 
cursor positioned over the ''window in'! command. (From Ramsden 


et abe: 1976). 
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FIDUCIALS OFF ON 
WINDOWS | OFF ON 
LANDMARKS OFF ON 
CONTOURS OFF ON 
GRND REF SET CLR 
OUTPUT OPTIONS: 
HARD COPY OFF ON 


WINDOW OFF ON 
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COMMENT 
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is able to select a portion of the PPI (a "window") to extract data for numerical 
processing. Figures 4.2 and 4.3 depict the "window" feature. Subsequent dis- 
plays will consist only of those data lying within the indicated portion of the PPI. 

The radar echoes can be represented by radial lines or by contours of 
equivalent radar reflectivity using the "contours on" display option (Figure 4.4). 
The lowest contour value is always 20 dBZ; the contour interval may be chosen by 
the analyst to be any value from 0.5 dBZ to 20 dBZ. 

Using this technique, individual storms were tracked and the data 
extracted from which the values for 7 were calculated. 

4.2 Storm Selection 

The analysis was aimed at determining what effects cloud seeding 
may have on the parameter 1 and if indeed the parameter 1 is correlated with 
hailfall at the surface to a degree sufficient to warrant its use as an evaluation 
parameter. Therefore it was necessary to choose storms where the time of seeding 
and the extent of hailfall were accurately known. 

Only storms with well-defined boundaries during most of their lifetime 
were chosen. In general, a storm was chosen if a window could be drawn around 
the storm echo, and tracked in this way throughout its' lifetime within the project 
area, with no other storms intruding or merging with it. This eliminated storms 
which were parts of "squall lines" or of agglomerations of many embedded storms. 
All storms selected were from the northern "randomized" area of the project. Only 
those which tracked sufficiently distant (more than 36 km) from the radar antenna 


thus permitting near constant viewing aspect were used. Both seeded and non-seeded 
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storms were chosen in this manner. 


4.3 Data Sample 


For the purposes here, 10 storms were chosen for analysis. Of these, 


5 storms were seeded and 5 storms were not seeded. Hail was reported to have 


fallen from 9 of the storms, with only rain reported from the remaining storm. The 


storms will be referred to by date, with letters A, B, etc. added to distinguish 


several storms occurring on the same day. Table 4.1 lists the entire storm sample 


for the analysis. 


Table 4.1 Storm Sample 


No. 

ig 20 JULY 1975: A 
Zz 20 JULY 1975: B 
eg 20 JULY 1975: C 
4, 17 JULY 1976: A 
53 17 JULY 1976: B 
é. 6 AUGUST 1976 
ve 12 AUGUST 1976 
8. 15 AUGUST 1976 
9. 20 AUGUST 1976 
10. 24 AUGUST 1976 


Rain: 
Hail: 
Hail: 
Hail: 
Hail: 
Hail: 
Hail: 
Hail: 
Hail: 


Hail: 


non-seeded 


seeded 


seeded 


seeded 


non-seeded* 


seeded 


non-seeded 


non-seeded 


seeded 


non-seeded 


* 1 cloud top flare was dropped into the storm early in its lifetime but for this 


analysis it has been considered to be non-seeded. 
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The 7 storms from the 1976 hail season were the only storms in 1976 
which met the requirements for selection, therefore any subjective bias in sample 
selection in these cases is minimal. 

The 3 storms of 20 July 1975 were chosen mostly for comparison 
purposes. As will be discussed later, the 1975 digitized radar data was often of 
questionable quality with certain portions of it missing. 20 July, 1975 was a rare 
occasion in that 3 storms met the requirements for selection thus, the storm data 
could be compared with other data on the same day as well as the 1976 data, to 
see if the quality was consistent and acceptable for this type of analysis. 

The data were extracted and the values of " were calculated. In 
most cases a PPI was obtained every 1.5 minutes. However, certain portions of 
the storms on 12 August 1976, 15 August 1976, and 24 August 1976, occurred while 
the radar was on a full 20° elevation program and, therefore, a PPI at the desired 
level occurred at 3.0 min intervals. 

4.4 The Relation Between 1 and Hailfall at the Surface 

In order to use the parameter 7 to test for cloud seeding effects, it 
was first necessary to verify its' assumed relation with hailfall at the surface. Res- 
ponses from farmers using hail reporting cards, and information solicited through 
telephone surveys following storms, were the major sources of hailfall data. 

Since radar reflectivity (and consequently the evaluation parameter) 
vary as the sixth power of the target diameter and only as the first power of the 
number concentration, it was decided to relate the reported maximum size of the 


hail which fell at the surface to the value of which existed at that time. Although 
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the time of the echo for which 7 is calculated is known to within an accuracy 

of several seconds, the time of a hail report by a farmer is unlikely to have a 
precision better than 5 minutes (Williams and Douglas, 1963). Moreover, because 
certain regions within the project area, particularly the western portion, are 
sparsely populated, hail reports may not have been received from regions where hail 
fell. Furthermore, even if a hail report is correct, the comparison with 7 is 
accurate only if the point surface observation is representative of the precipita- 
tion which fell from the volume sampled by the radar rn: 

The method used to test the relation of 1 with hailfall at the surface 
was to combine the maximum hailsize reports with the values of 1 recorded 
closest in time for all 9 hailstorms. The time involved for the hail to fall from 
the echo volume to the ines was not accounted for in view of the accuracy of 
the reports. Although some reports are in error and contribute to the variance, by 
combining all of the reports and assuming more good reports than bad reports, the 
averages are likely to indicate the true relation. 

The farmers are asked to report the maximum hailsize in one of 6 
categories associated with the sizes of familiar objects (Table 4.2). A total of 358 
hail reports were received from the 9 hailstorms that were analysed. Figure 4.5 


shows the cumulative distribution of all the hail reports according to their 1 values. 


The range of 1 values for which hail was reported was from 0.62 x 10-!!em-!to 
2.57 x 10? cm~} with a median value of 4.20 x 101% cm!. 
Figure 4.6 shows the cumulative frequency curve for the individual 


maximum-size categories. The categories for golfball and larger-size hail were 
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Table 4.2: Hailsize name and the assumed hail diameter range. 


Hailsize Name Assumed Diameter Range (cm) 
Shot “i 6.3 
Pea 4— 1.2 
Grape rod 2.0 
Walnut Palio. 
Golfball = ie ot ae 
Larger as PY 


Source: Wojtiw, 1975b 


combined since only | report of hail larger than golfball size was reported. One 
can see that, within each size category, a large spread in 1 values occurred. 
The larger 1 values within the smaller hailsize categories may be due to the 
fact that the farmer may be located at the edge of the hailswath and receiving 
small hail even though the storm may be producing large hail which is falling 
elsewhere. However, if 7 is correlated well with hailfall, large hail should 
not be associated with small values of 7. Walnut size hail was reported at quite 
low values of 7, but these reports are few and may well be ones where the time 
is in considerable error, or the maximum hailsize was overestimated. 

A plot of the mean 7 values for the 5 categories of maximum hailsize, 
along with the 90% confidence intervals for the individual means calculated using 


the statistic t-test, is given in Figure 4.7. The median and mean 7 values for 


the individual size categories indicate that larger hail tends to occur with increasing 
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Fig. 4.7: The mean n values obtained for the five categories of maximum 
hailsize, and the 90% confidence interval. 
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1 values for the maximum size categories shot through walnut. A deviation in 
the trend occurs for the hailsize category of golfball and larger. There isa 
sampling problem relating to the largest hail sizes since only 16 reports of golf- 
ball and larger size hail were received; however, this trend of increasing, then 
decreasing reflectivity with size has been found by others. Theoretical calcula- 
tions for monodisperse hail by Atlas and Ludlam (1961) showed that radar reflec- 
tivity has an oscillatory behavior once the size of the scattering particles extends 
beyond the Rayleigh region (Figure 4.8). The maxima and minima in reflectivity 
for scatterers which are no longer small with respect to the radar wavelength 

are characteristic of Mie's scattering functions. Sulakvelidze et al. (1965) ob- 
tained smooth curves for the diameter versus radar reflectivity relationship for a 
series of calculated Harlefene spectra (Figure 4.9). The N versus maximum hail- 
size relation summarized in Figure 4.7 resembles the curves of the theoretical 
work cited. Therefore, the decreasing values of 1 associated with hail larger 
then 3 cm diameter are not likely due to a sampling problem, but rather, due to 
the inherent limitations of radar to adequately differentiate between large scat- 


terers. 


It is worthy to note the merit of the 10 cm wavelength radar when 
viewing hail. Shorter wavelength radars yield oscillating radar reflectivities at 
smaller target diameters and also attenuate severely, thus making hail detection 
more difficult. 

The comparison of maximum hailsize to 7 values in this case using 


a 10.4 cm wavelength radar, indicates that 7 is related to hailfall at the surface 
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Radar reflectivity as a function of the diameter of uniform, 
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(2) 5m-3; (3) im-3; (4) 0.01m™2. (From Sulakvelidze et al.,1965). 
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such that the probability of larger size hail falling at the surface increases with 
increasing values of 7 for hailstone diameters less than 3 cm. 
4.5 Time Series Analysis 

For each storm the 7 values were plotted as a function of time. For 
seeded storms, the type, amount and time of seeding were also noted. Although 
the detailed fluctuations of a time series may be quite irregular, certain statis- 
tical properties may remain fixed from one period to another. A statistical analysis 
was performed on each time series, to better understand the basic properties of 
the time series, their variability, and the characteristics of their periodic and 
irregular oscillations. Since the time series variate 7 is related to the type, size 
and amount of precipitation, such an analysis helps in the physical understanding 
of storms. 

The time series analysis consisted of computing the power spectrum 
by the fast Fourier transform (FFT) method. The computer program used was one 
developed and tested by Lachapelle (1977). The FFT is a highly efficient procedure 
for computing the discrete Fourier transform (DFT) of a time series consisting of 
discrete data samples. Its advantage lies in the fact that the coefficients for the 
DFT are calculated iteratively, which results in a considerable saving of compu- 
tation time. Given a record of N observations taken at a set time interval At, 
a total of N/2 harmonics can be calculated, the last one having a period of 2At. 
The contribution of each harmonic to the total variance of the time series is refer- 
red to as the power of each harmonic, and the power spectrum of a time series shows 


the contributions of oscillations of different frequencies to the variance of the time series. 
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Prior to the FFT analysis, the program required that the data be 


modified accordingly: 


ic The mean M, variance S*, and standard deviation S$, were calculated. 

oe The mean was removed. 

3.  Zeroes were added to the data set until a desired size of 2" was 
achieved. 

4, The data were standardized using the relationship 


Vio (kes M)/S 

Behe b ite a? let ped se 

Several reasons for statistical spectroscopy on meteorological variates 
are given by Panofsky and Brier (1963): 

1. The spectrum may be necessary to understand the physics underlying 
the variation of a time series. 

2. Significant maxima and minima are important for forecasting, even 
if they are not sharp peaks or troughs, for they indicate the likeli- 
hood or unlikelihood of variations with certain average periods. 

3. Cloud seeders and other personnel engaged in weather control claim 
they can change the spectrum significantly. 

4. The extent of the spectrum shows the variability of the quantity and 
how quickly instruments have to respond in order to measure it. 


The results of this analysis follow. 
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CHAPTER 5 
THE 17 JULY 51976 GASE STUDIES 

5.1 Introduction 

Two severe hailstorms occurred on 17 July, 1976. Figure 5.1 depicts 
the synoptic situation which prevailed and includes a brief synoptic summary. 
Figure 5.2 shows the rawinsonde sounding from Calgary at 18:00 MDT which de- 
scribes the environmental conditions. The two storms formed in approximately the 
same region and tracked approximately across the same area. What is of particular 
interest is that one storm was heavily seeded while the other was essentially un- 
seeded (only one cloud-top flare released into it), because of the late hour at 
which it entered the project area. These two storms resulted in 242 reports of hail 
at the surface acceptable for this study. This represents 67.6 per cent of the total 
number of hail reports used in determining the relation of 7 with hailfall at the 
surface. Therefore, the relationships deduced in Chapter 4 were greatly influenced 
by what occurred on this day. 

The radar constant used for the day was 73.2 dBm with a peak trans- 
mitted power of 286 kw. | 
Oe2 Ilo) 276 Storm A 

Storm A entered the northwest portion of the project area shortly after 
19:00 MDT and tracked eastward. Figure 5.3 shows the resulting hailswath of 
storm A. Hail as large as golfballs was reported at the surface. The dashed line 
attempts to isolate the region of hail from the region where rain only was reported 
to have fallen. The approximate times of storm progression are also indicated. 


Radar-estimated storm tops exceeded 10 km AGL, with radar reflectivities greater 
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SYNOPTIC SUMMARY 


17 JULY 1976 


SYNOPSIS: AFTER 14 JuLy, L/W PATTERNS SHIFTED W TEMPORARILY. WITH THE 
L/W RIDGE OVER ALBERTA FOR A FEW DAYS, S/W'S WERE DRIVEN NORTH OF ALBER- 
TA AND A PRONOUNCED SURFACE HIGH PRESSURE CELL BUILT OVER WRN CANADA, 
THIS L/W RIDGE THEN MOVED E ONCE AGAIN ALLOWING THE FIRST OF THE NEW 
SERIES OF S/W TROUGHS TO REACH THE AHP AREA LATE ON THE 17TH. ALTHOUGH 
NOT A WELL MARKED TROUGH (NEGLIBLE PVA FLOWS), 500 MB HEIGHTS FELL RA- 
PIDLY. MEANWHILE THE SURFACE HIGH DROPPED SE AND A WAVE FORMED ON THE 
FRONT IN SE B.C. THE COMBINED EFFECTS TRIGGERED THE WAITING INSTABILITY 
AND ONE OF THE LARGEST STORMS OF THE YEAR ENSUED, 


AHP AREA: LIGHT NE SURFACE WINDS VEERING TO STRONG SW AT 500 mB. Two 
LARGE STORMS FORMED RAPIDLY IN THE NORTH DURING EARLY EVENING, PROPA- 
GATING ESE WITH TOPS EXCEEDING 12 kM, 


12 HOUR 500 MB CHANGE IN 


VORTICITY|HEIGHTS *{THICKNESS*| TEMPS. 
(o> °s: 1) | tu) (M) (°C) 


TROUGH (T) OR 
RIDGE (R) 
PASSAGE TIME 
R-17/0100Z 


T-18/06002 


HAIL 
SLZE 
CATEGORY 


NO. HAIL 
REPORTS 
Naes 


S00 MB ANALYSIS y 
HEIGHTS AND VORTICITY 
0000Z 18 JULY 1976 


Fig. 5.1: Map showing the synoptic situation, and a brief synoptic 


hy 


summary for storms on 17 July, 1976. (From Deibert (ed.), 1977). 
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STATION: CYYC 
TEPHIGRAM DATE: 18 JULY 76 


TIME: 00002 


Rawinsonde sounding from Calgary at 18:00 MDT 17 July, 1976. 
(From Deibert (ed.), 1977). 
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than 60 dBZ. Figures 5.4 through 5.8 show low-elevation PPI scans at 
approximately 15 minute intervals throughout the storm lifetime. 

Figure 5.9 shows the 7] time series as computed from the digitized 
radar data. The right-hand vertical axis gives the equivalent radar reflectivity 
values computed from 1 using (2.2.4). In calculating Za, the quantity [e+ 2/e-1]* 
conventionally is taken to be the value for water 0.93 (Battan, 1973). The histo- 
gram at the top of the graph indicates the number of cloud-seeding flares released 
during each 15 minute interval as reported in the flight operations summary for the 
day. The dashed line superimposed on the time series is the trend in the data 
determined by using a low-pass Tukey filter described in Appendix A. 

It is of particular interest to note that the n value at 21:00 MDT is 


“1 arelatively low value, whereas the PPI scan at 


approximately 1.6 x 107° cm 
21:00:41 MDT indicates that reflectivity exceeded 60 dBZ. Although one would 
not expect as low an 7 value with such high reflectivity, it is apparent from the 
PPI that the reason is that high reflectivity volumes have decreased while the over- 
all storm extent, bounded by the 20 dBZ contour, has increased slightly. The time 
21:00 MDT corresponded with the break in the hailswath at the surface which in- 
dicates that the storm intensity decreased significantly at this time. This would not 
be deduced readily from the PP! alone. In fact, one may be inclined to suspect the 
storm to have intensified since the areal extent of the echo had increased. In 


this instance the parameter " appears to be quite sensitive to storm intensity. 


5.3 Storm A: Power Spectrum 


The power spectrum of the raw data is presented in Figure 5.10. Although 
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Fig. 5.4: PPI echoes for storm A on 17 July, 1976. Elevation ~ 250. 
Reflectivity = from 20dBZ by 10dBZ. Time = 19:30 - 20:15 MDT. 
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Fig. 5.5: PPI echoes for storm A on 17 July, 1976. Elevation ~ a 
Reflectivity = from 20dBZ by 10dBZ. Time = 20:30 - 21:00 MDT. 
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PPI echoes for storm A on 17 July, 1976. 
Reflectivity = from 20dBZ by 10 dBZ. 
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hia. Boas: PPRPI «echoes for sstormfAvon (17 July, 1976. iGlevatiion ~~ Dor 
Reflectivity = from 20dBZ by 10dBZ. Time = 21:45 - 22:00 MDT. 
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Fig. 5.8: PPI echoes for storm A on 17 July, 1976. Elevation ~ 2.2° 
Reflectivity = from 20 dBZ by 10 dBZ. Time 22:15 - 22:45 MDT. 
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8 ee tins J n versus time for storm A on 17 July, 1976. Dashed line indicates 
trend. Type and number of cloud-seeding flares released during 
15 minute intervals is shown. 


50 


oe 


oH sibal ane cain a Sie ee re ne Wa 
eh ul i al ne ral wa obits, 0 


| - . 4a os | oi) a 
} yy, j } - 7 ae , : 
ele: a a oat te © as 


51 


2 
1.076 
0.000 0.033 0.066 0.099 0.132 0.165 


CYCLES PER MINUTE 


Fig. 5.10: Unsmoothed power spectrum of n for storm A on 17 July, 1976, 


showing the white-noise spectrum and 90% confidence interval. 
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Fig. 5.11: Three-point Daniell-filtered power spectrum of n for storm A 


on 17 July, 1976, showing the red-noise spectrum. 
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all possible harmonics were computed, (up to 0.333 cycles per minute), the spec- 
trum is only plotted for frequencies up to 0.165 cycles per minute or periods of 
approximately 6 minutes (4 At). Figure 5.11 shows the power spectrum which has 
been smoothed using a three-point running average applied to the periodogram. 
The reason for smoothing is that, if the spectra of two different portion of the same 
stationary time series are derived, the individual points of the spectrum will most 
likely look quite different. Spectrum analysis does not attempt to determine the 
individual amplitudes of each harmonic, but is aimed at obtaining the smooth under- 
lying spectrum which is the same for different portions of the same stationary time 
series. It was once thought that particular outstanding powers represented exact 
periodicities, however, it has been found that individual peaks in the spectra of 
irregular oscillating time series can be accidental at a certain frequency. 

The three-point running average also called the Daniell power spectral 
estimate (Kanasewich, 1975), produces an estimate of the continuous spectrum 
viewed through a rectangular window of bandwidth equal to 3 f/m, where f,, is 
the folding or Nyguist frequency (fy = 1/2 At), and m is the number of spectral lines. 
The end points were ignored with the filter applied only os needed. Many authors 
(Blackman and Tukey, 1958; Jenkins, 1961) express the need for bringing the re- 
sultant spectrum close to that of "white noise"; also called "prewhitening" or the 
removal of the trend in the data. For a full explanation of white noise see Appendix 


B. The reason for prewhitening is that the average value of the power at any 


particular frequency may be distorted during computation since the effect of a spectral 


window is to spread the power from the large peaks into adjacent frequencies. Thus 
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spectral estimates are considered to be more precise when the power is spread over 
all frequencies. 

The power spectrum was computed using the Tukey filtered high-pass 
data and is shown in Figure 5.12. It was found that the FFT program was able to 
discern the high frequency powers capably in the presence of substantial low 
frequency power, nevertheless, a smoothed FFT spectral analysis of the high pass 
data did bring out more clearly any preferred small scales. 

5.4 Storm A: Spectrum Analysis 

The value of the white-noise power is applied to the unsmoothed power 

spectrum in Figure 5.10 assuming a Chi-square distribution. The number of degrees 


of freedom v for the periodogram analysis was determined by using: 


ie eee (52451) 
ie 


where N? is the number of independent data values, 2°~ * is the number of spectral 


lines obtained from N data values supplemented by 2" - N zeroes, and D is the 
width of the Daniell filter used on the periodogram. The number of independent 
data values in the original N pieces of data depends on the amount of autocorrela- 


tion ©, and has been shown by Eddy et al., (1968) to equal to: 


1- 0? 
UD Peet 
Di y(24) 


The number of independent data values is then reduced by 2 degrees of freedom since 


the mean and variance were calculated from the data. 


noo 


i 


ct ena ane Sobge yaw 
+ pee. 
et siesta en fh soa = 


‘ 
‘1 * 


. wat itt ig ln or 
7 us ons he nee hh 


Mee = ba diccagru aft a lle a wuewreny scones ont Va ve | 
ob Ty nechown it -eai' 7 hte aninepa a: satoonsDE2 at 
si teuiad Geel mr oiit Boa adlit ec chen a 


ait 


in $ 3 
j 7 _ si _ “4 te 
iooge ta stdin edt a 2 Qteuliay oe roabrernbed Yo rdkewn at ah 


\ 
‘ 


we at tl wore es Me "Sb vit ai sale sah Hc 


Poy 


at = a - 4 
verb agabrsi te ebaiaall eii| Mey rrr we, ne be sai ae . hey: 
or in he” Sas rae: . 


olavenss oS aes te deniomiy ext. no surmgd ooh te sonata v izoighe a 
ot oe ro (Get) ae 0 wn rel ect 


aunts @nbeatt ‘ere Ft 


| | | ‘ee rs! eo eperre 
m i i) 7. - a _ - i Proceed ) 
a peo 


a 


} | f : 7 ‘“ - : if : 
hot nen 


POWER 


PUG ow Zc 


0.000 0.033 0.066 0.099 0.132 0.165 
CYCLES PER MINUTE 


Three-point Daniell-filtered power spectrum of Tukey-filtered 
n highpass data for storm A on 17 July, 1976. The red-noise 
spectrum and 90% confidence interval are shown. 
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The 95% confidence limits were obtained by multiplying the white- 
noise value by x* .95/v and x*.05/ v for the upper and lower limits, respectively. 
The 1 series does not display random fluctuations characteristic of white noise. 

The spectrum displays a general suppression of relative variance at 
higher frequencies and consequent inflation at lower frequencies as compared to 
the even distribution of white noise. This phenomenon, a characteristic of spectra 
derived from data with an appreciable autocorrelation between successive mea- 
surements, has been called "red noise" (Gilman et al., 1963). For a full explana- 
tion of red noise see Appendix B. 

The lag one (1.5 min) to lag ten (15 min) autocorrelation coefficients 


were computed as follows: 


N- 7 
ROT) pee eee 2X CX) MOG ae =X) (5.4.2) 
IN Ay Se See 


where N is the number of data points (reduced by 2 in denominator since the mean 
X, and the variance SS were computed from the data), and T is the lag number. 
Table 5.1 gives the results. 

The autocorrelation coefficients confirm that the time series depends 
on its immediate past and that a red-noise spectrum which allows for persistence 
should be applied. 

The values of the red~noise power are applied to the smoothed power 
spectrum in Figure 5.11 using an autocorrelation value © = .992 for the red-noise 


curve determined from the relation po = [R(1)® + R(2) }3/7 . From the definition 
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Table 5.1: Lag 1.5 to 15 minute autocorrelation coefficients of n 
for~stormtAeonelginulyeril976: 


LAG NUMBER LAG TIME (min) AUTOCORRELATION COEFFICIENT 


1.5 999 
2 3.0 970 
3 4.5 936 
r 6.0 898 
5 7.5 855 
é 9.0 .810 
7 | 10.5 765 
8 1220 719 
9 ial es .672 
10 15.0 625 


of pure red=noise 


Since this is not the case with experimental data, the value of 9 should be mod- 
ified to approximate pure red-noise. It was this value for © which was also used 
in the calculation of N°. 

The values of the red-noise power are applied to the smoothed high- 


pass power spectrum in Figure 5.12 using an autocorrelation value 9 = .684. We 


may assume a Chi-square distribution with the number of degrees of freedom v 


calculated using (5.4.1). In this case D=3, since a three-point running average is used. 
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The red-noise adequately describes the nature of the smoothed power 
spectrum. The computed spectral estimate does not display any significant peaks 
with respect to the 95 per cent upper confidence limit, nevertheless, preferred 
cycles were displayed having periods of 32.0, 16.0 and 9.6 minutes. 

The possible relationships between the variation of the time series and 

cloud-seeding events will be discussed in a separate section after all the indi- 
vidual storm cases have been presented. 
5.5 17 July, 1976: Storm B 

Storm B entered the project area shortly after 21:00 MDT and tracked 
eastward just south of storm A. It also had echo-top heights similar to those of 
storm A. Figure 5.13 shows the resulting hailswath from storm B with approximate 
times of the storm progression indicated. Golfball size hail and one report of 
larger-than-golfball size hail was reported at the surface. Figure 5.14 through 


5.18 show low-elevation PPI scans at approximately 15 minute intervals throughout 


the lifetime of storm B. Once again, reflectivities in excess of 60 dBZ were measured. 


Figure 5.19 gives the 7 time series. Only one cloud-top seeding 
flare was released into the storm at 22:05 MDT, as shown in the diagram. The dashed 
line indicates the trend in the data computed using the Tukey filter in Appendix A. 
A notable feature of the time series is that it displays a relatively uniform nature 
as compared with storm A which peaked, decreased significantly in magnitude and 
then surged to a peak once again before dissipating. Storm B reached a high plateau 
which gradually peaked at approximately 23:30 MDT and then dissipated. This 


difference is also displayed in the form of the hailswath. The hailswath for storm A 
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PPll echoes’ for storm open lv July, 
Reflectivity = from 20dBZ by 10dBZ. 


1976. 
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Fig. 5.15: PPI echoes for storm B on 17 July, 1976. Elevation ~ one 
Reflectivity = from 20dBZ by 10dBZ. Time = 21:45 - 22:30 MDT. 
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Fig. 5.16: PPI echoes for stornmB on 17 July, 1976. Elevation ~ 2.] 
Reflectivity = from 20dBZ by 10dBZ. Time 22:45 - 23:30 MDT. 
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Fig. 5.17: PPI echoes for storm B on 17 July, 1976. Elevation ~ 2.2° 


Reflectivity = from 20dBZ by 10dBZ. Time = 23:45 - 00:29 MDT. 
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Figw oerocs Pri, echoes for stotmecson 17 July, 1976> Elevation ~ ORs 
Reflectivity = from 20dBZ by 10dBZ. Time = 00:44 - 01:29 HoT. 
(18 July, 1976) 
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Fig. 5.19: n versus time for storm B on 17 July, 1976. Dashed line indicates 
trend. The time when one cloud-top flare was released in the 
storm is noted. 
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displays an intermittency while the hailswath for storm B is continuous. The 7 
values obtained for storm A at its onset were considerably higher than any values 
obtained for storm B. Although this difference is not evident in the size of hail 
reported at the surface early in storm A, this is due possibly to inadequate sampling 
in a sparsely-populated area. Comparable precipitation did occur otherwise for 
similar 7 values. Storm B continued on into the early morning hours of the fol- 
lowing day. 

5.6 Storm B: Power Spectrum Analysis 

Figure 5.20 gives the unsmoothed power spectrum for storm B. The 
smoothed power spectrum using a three-point running average is presented in 
Figure 5.21. 

As was the case with storm A, the power spectrum displays the red- 
noise phenomenon. Table 5.2 gives the autocorrelation coefficients of the time 
series for lag one (1.5 min) to lag ten (15 min), computed using (5.4.2). 

The values of the red-noise power are applied to the smoothed power 
spectrum in Figure 5.21 using an autocorrelation value 9 = .986, and to the 
highpass data smoothed power spectrum in Figure 5.22 using an autocorrelation 
value ©9= .299. The 95 per cent confidence limits were constructed assuming 
a Chi-square distribution as in storm A. 

Significant peaks exist at 0.049 and 0.109 cycles per minute cor- 
responding to cycle periods of 20.4 and 9.2 minutes respectively. 

5.7 Spectra Comparison 


Storm B displayed a noticeable lack of power at 0.066 cycles per minute 
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Fig. 5.20: Unsmoothed power spectrum of meforestorm 6 only July. 19/6. 
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Fig. 5.21: Three-point Daniell-filtered power spectrum of n for storm B 
on 17 July, 1976, showing the red-noise spectrum 
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Three-point Daniell-filtered power spectrum of Tukey-filtered 
n highpass data for storm B on 17 July, 1976. The red-noise 
spectrum and 90% confidence interval are shown. 
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Table 5.2: Lag 1.5 to 15 minute autocorrelation coefficients of 1 
for storm Blon-d7eJuly, 1976: 


LA G: NUMBER TIME LAG (min) AUTOCORRELATION COEFFICIENT 


ies 990 
2 3.0 963 
3 At5 94] 
4, 6.0 926 
5 7.5 913 
6 9.0 893 
w) 7 10.5 .866 
8 1220 .834 
9 ele-5 .805 
10 15.0 .783 


whereas storm A displayed a preferred frequency. Storm A does, however, show 
a broad peak coinciding with the significant peak for storm B at approximately 
0.11 cycles per minute, although well within the confidence interval. 

The Chi-square test for goodness-of-fit was applied to the two raw 
spectra to see if the observed power distribution for storm A differed significantly 
from that of storm B. The results indicated that even at very low confidence limits ~ 
the two spectra could constitute a sample from a population with the same power 
distribution function. 

A physical interpretation of the preferred cycles exhibited is discussed 


in a later section after all the storm data has been presented. 
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CHAPTER 6 
THE AUGUST 1976°GASE STUDIES 

6.1 Introduction 

Five storms which occurred during the month of August, 1976 were 
analyzed. The storms chosen on August 6, 12, 15, 20 and 24, combined to yield 
67 reports of hail at the surface and are included in this study in order to compare 
the variations of 7 in storms with lower energy than the energy associated with 
the two storms discussed in Chapter 5. 

The storms on August 6 and 20 were seeded while those on August 12, 
15 and 24 were not. Since smaller storms are generally thought to be among the 
most likely to be affected by cloud seeding, seeding effects may be more easily 
detected. 
6.2 The 6 August, 1976 Storm 

The storm on 6 August, 1976 differs from most of the other storms in 
this study because it did not form in the foothill regions and track eastward across 
the project area, but rather formed on the plains in the vicinity of Buffalo Lake 
and went through its life cycle within a region of less than 300 km® over a time 
period of 1.5 hours. The synoptic summary is given in Figure 6.1 and the environ- 
mental conditions are shown by the tephigram in Figure 6.2, Although the presence 
of the lake would offer an opportunity to study the thermodynamic influence of a 
lake on a convective storm for comparison with the lake effect proposed by Warner 
(1976), this analysis was not undertaken. Before any changes in 1n can be associated 
with cloud seeding, however, the presence of the lake would have to be taken into 


consideration. 
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SYNOPTIC SUMMARY 


06 AUGUST 1976 


SYNOPSIS: THE SECOND S/W TROUGH IN THE CURRENT SERIES ARRIVED OVER THE 
AHP AREA AROUND NOON, THE SURFACE PRESSURE PATTERN DID NOT REFLECT THE 
S/W TOO WELL,WITH ONLY A WEAK FRONTAL TROUGH ACROSS CENTRAL ALBERTA, 


AHP AREA: LIGHT AND VARIABLE WINDS AT ALL LEVELS. MOST CONVECTIVE DE- 
VELOPMENT OVER THE AHP AREA TOOK PLACE IN THE NE QUADRANT, WHERE THE 
TIMING OF THE SYNOPTIC S/W SUPPORT WAS MOST FAVORABLE. CELLS FORMED 
RAPIDLY AROUND NOON, WITH SOME TOPS EXCEEDING 13 KM BY MAXIMUM TEMPERA- 
TURE TIME WHILE MOVING NORTHWARD. 
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Fig. 6.1: Map showing the synoptic situation and a brie 
for phavetormscn-6 August 91976+ ‘(From Detbert (ed.), 1976) 
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7 August, 1976. (From Deibert (ed.), 1977) 
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The hailswath from the storm is shown in Figure 6.3. Only 5 reports 
of hail at the surface were received and the maximum hailsize reported was grape 
size. 

PPI scans of the storm at approximately 15 minute intervals are given 
in Figure 6.4. The storm developed from the minimum detectable signal to reflec- 
tivities greater than 50 dBZ in less than half an hour. 

The 7 time series is given in Figure 6.5, which also indicates the 
time and amount of seeding which was done. The dashed line represents the trend 
in the data computed using the Tukey filter described in Appendix A. 

The radar constant used for the day was 75.5 dBm with a peak trans- 
mitted power of 168 kw . 

The unsmoothed power spectrum of the data is given in Figure 6.6 and 
the smoothed power spectrum, using a three-point running average, is given in 
Figure 6.7. The frequency resolution of the spectrum is not as good in this case 
as in the previous cases. The storm lasted just over 1.5 hours and the number of 
data points available was 63. The FFT algorithm requires that the number of data 
points be a power of two. If the number of values is not a power of two, zeroes 
are added until the total number meets this requirement. The number of spectral 
lines between frequencies 0.0 and (1/2 At) cycles per minute depends on the rela- 
tion m = 2°-+ where m is the number of spectral lines and n is an integer. In this 
case the value for n is 6 and therefore, Vee 64, which means that one zero was 
added to complete the data set. The number of spectral lines in this case is 32. 


The number of spectral lines for each storm on 17 July 1976, was 128 resulting in 
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Fig. 6.3: Map showing the hailswath resulting from the storm on 6 August, 1976. 
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Fig. 6.4: PPI echoes for the storm on 6 August, 1976. Elevation ~ yikes 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 16:31- 17:59 MDT. 
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indicates trend. Type and number of cloud-seeding flares released 
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Unsmoothed power spectrum of n for the storm on 6 August, 1976. 
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Three-point Daniell-filtered power spectrum of n for the 
storm on 6 August, 1976, showing the red-noise spectrum and 
90% confidence interval. 
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considerably more detail. This is the reason for the very smooth looking spectrum 
of 6 August, 1976 as compared with those from 17 July, 1976. 
Table 6.1 gives the autocorrelation coefficients of the time series for 


lag one (1.5 min) to lag ten (15 min), computed using (5.4.2). 


Table 6.1: Lag 1.5 to 15 minute autocorrelation coefficients of n 
for, the stonm.on 6 August, 1976. 


LAG NUMBER TIME LAG (min) AUTOCORRELATION COEFFICIENT 


1 ae ae 
a | 3.0 891 
3 4.5 a7 07 
4, 6.0 654 
5 4 so27 
6 7.0 430 
i LOe5 O20 
8 12.0 238 
9 1325 .140 
10 15.0 054 


The values of the red-noise power assuming a Chi-square distribution 
are applied to the smoothed power spectrum in Figure 6.7 using an autocorrelation 
value P= .958. 

There are no small-scale preferred frequencies consistent throughout 
the data. The application of the Tukey highpass filter used in this study results 


in 10 data points being lost at each end of the data set. This loss is unacceptable 
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in this case due to the small sample size originally, therefore the power spectrum 
of the highpass filtered data revealed no useful information. 

It should be hx that the two cloud=base flares were burned at 17:22 
MDT and 17:27 MDT, respectively, and that 8 cloud-top flares were released at 
17:29 MDT according to the flight summary. The decline in 7 preceded the 
cloud seeding. This is not obvious from Figure 6.5 which does not show the exact 
time of seeding. 

The storm on 6 August, 1976 exemplifies the sampling problem involved 
in using only surface hailfall measurements in evaluating the output of a given 
storm since the presence of the lake for this case severly limited the number of 
observations obtainable. 

6.3 The 12 August, 1976 Storm 

The storm on 12 August, 1976 entered the project area from the north- 
west and spent over five hours in the northwest quadrant of the project area before 
dissipating. The storm was not seeded. 

The synoptic summary is given in Figure 6.8 and the environmental 
conditions are shown by the tephigram in Figure 6.9. The hailswath is shown in 
Figure 6.10 which indicates that 6hailreports were received, including one report 
of golfball size hail. It is difficult to say if this report was accurate since so few 
hail reports were received from this sparsely-populated area. 

PPI scans of the storm at approximately 15 minute intervals are given 
in Figures 6.11 through 6.14. Radar reflectivities exceeded 40 dBZ at the time 


the storm entered the project area. The storm was very slow moving. Two centers 
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SYNOPTIC SUMMARY 


12 AUGUST 1976 


SYNOPSIS: THE L/W RIDGE NOW SHIFTED E OVER SASKATCHEWAN. A N-S SURFACE 
TROUGH OF LOW PRESSURE FORMED OVER ALBERTA IN RESPONSE TO A S/W TROUGH 
PASSING JUST NW OF THE AREA; THIS LED TO AN INCREASE IN LOW-LEVEL MOIS- 
TURE OVER THE PROJECT AREA, 


AHP AREA: LIGHT EASTERLY SURFACE WINDS BECOMING SW AT 500 MB. THE MOST 


SIGNIFICANT STORM CELL DEVELOPMENTS WERE ALONG THE FOOTHILLS WHERE MAXI- 
MUM TOPS EXCEEDED 13 KM BY LATE EVENING WHILE MOVING N, 


TROUGH (T) OR 12 HOUR 500 MB CHANGE IN HAIL |NO. HAIL 
RIDGE (R) VORTICITY |HEIGHTS *|THICKNESS*| TEMPS.| SIZE REPORTS 
PASSAGE TIME (Tepes) (M) (M) (OMe PCATEGORY |) aN fas 


T-13/01002Z 
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: /\ BS 4 
500 MB ANALYSIS 
HEIGHTS AND VORTICITY 
00002 13 AUG 1976 


Fig. 6.8: Map showing the synoptic situation, and a brief synoptic 
summary for the storm on 12 August, 1976. (From Deibert 


(ed.), 1977) 
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Fig. 6.9: 


fey : 
O STATION: CYQF 
TEPHIGRAM DATE: 12 AUG 76 


TIME: 21002 


Rawinsonde sounding from Red Deer at 15:00 MDT, 12 August, 1976. 
(From Deibert (ed.), 1977) 
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Fig. 6.10: Map showing the hailswath resulting from the storm on 12 August, 1976. 
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Fig. 6.11: PPI echoes for the storm on 12 August, 1976. 


Elevation ~1.9° 
Reflectivity = from 20 dBZ by 10 dBZ. 


Time = 17:00 - 18:14 MDT. 
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Fig. 6.12: PPL echoes for the storm on 12 August, 1976. Elevation ~ 2.0° 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 18:29 - 19:45 MDT. 
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Fig. 6.13: PPI echoes for the storm on 12 August, 1976. Elevation ~ 2.0° 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 20:00 - 20:45 MDT. 
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Fig. 6.14: PPI echoes for the storm on 12 August, 1976. Elevation ~ 2.0° 
Reflectivity = from 20 dBZ by 10 dBZ. “Time = 21:00 - 22:00: MDT. 
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of precipitation are evident from the PPI scans over the period from around 18:00 
MDT to shortly before 20:00 MDT. After uniting, maximum radar reflectivity ex- 
ceeded 60 dBZ, and reflectivities exceeding 50 dBZ persisted until approximately 
21:30 MDT after which the storm dissipated rapidly. The 1 time series for 12 
August, 1976 is given in Figure 6.15. The dashed line represents the trend in the 
data computed using the Tukey filter described in Appendix A. 

The radar constant for the day was 74.2 dBm with a peak transmitted 
power of 227 kw. 

There existed during this day a period of time during which the radar 
was on a full 20 degree elevation scan. This meant that a PPI at one-degree ele- 
vation was obtained every 3 minutes instead of every 1.5 minutes as in the previous 
cases. The calculation of the power spectrum demands that the time interval remain 
constant between data points. Therefore, the data had to be modified in order to 
meet the requirement. Two alternatives existed. Firstly, every other data point 
during the 1.5-minute scans could be removed resulting ina data interval of 
3 minutes. Secondly, data could be interpolated between known data points during 
the time the radar was on the full 3-minute scan. The radar was put on the 3-minute 
scan at 21:00 MDT during this day. Since the storm ended at 22:03 MDT, only 
21 data points were missing. Therefore, it was decided that linear interpolation 
would be used. The error involved should be minimal since the storms already dis 
cussed had autocorrelation coefficients greater than .95 for a 1.5-minute time lag. 
Linear interpolation should represent a respectable approximation. 


The autocorrelation coefficients of the time series for lag one (1.5 min) 
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to lag ten (15 min) computed using (5.4.2) and including the 21 interpolated data 


points, are given in Table 6.2. 


Table 6.2: 


LAG NUMBER 


I 


2 


3 


4 


10 


Lag 1.5 to 15 minute autocorrelation coefficients of n 
for the storm on 12 August, 1976. 


TIME LAG (min) 


1 


3. 


4, 


6. 


ie 


Me, 


0 


§ 


0 


AUTOCORRELATION COEFFICIENT 
936 
821 
.687 
.554 
~457 
396 
365 
349 
.329 


31] 


The unsmoothed power spectrum of the data is given in Figure 6.16 


and the smoothed power spectrum, using a three-point running average, is given 


in Figure 6.17. The highpass data smoothed power spectrum is given in Figure 6.18. 


The values of the red-noise power are applied to the smoothed raw data power 


spectrum using an autocorrelation value 9 = .921, and to the highpass data smoothed 


power spectrum using an autocorrelation value 9 = .690. A broad peak at 0.057 


cycles per minute corresponding to a cycle period of 17.5 minutes appears signif- 


icant at the 95 per cent confidence level, assuming a Chi-square distribution in 
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Unsmoothed power spectrum of n for the storm on 12 August, 1976. 


O-000) nc.033 0.066 0.099 Hi0Gta2 0.165 
CYCLES PER MINUTE 


Three-point Daniell-filtered power spectrum of n for the storm 
on 12 August, 1976. The red-noise spectrum is shown. 
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Three-point Daniell-filtered power spectrum of Tukey-filtered 
n highpass data for the storm on 12 August, 1976. The red- 
noise spectrum and 90% confidence interval are shown. 
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Figure 6.18. It should be noted that at the 95 per cent confidence level, at least 
3 of the spectral lines in the 64 lines displayed in Figure 6.18 could reach statis- 
tical significance by chance due fo statistical fluctuations. Therefore, the broad 
peak in the power spectrum at 0.057 cycles per minute should more appropriately 
be interpreted as a preferred cycle. 

The time series displays greater small-scale variation and less auto 
correlation than the series previously discussed. This is most likely due to the fact 
that the storm was not a well-organized system but consisted of at least two preci- 
pitation centers, each with its own associated scales of variability, for a major 
portion of its lifetime. For example, fluctuations in the precipitation in one core 
would not coincide necessarily with the fluctuations in the other core. For this 
reason, this storm may not be a good example of a non-seeded storm for comparison 
with seeded storms. 

6.4 The 15 August, 1976 Storm 

The storm on 15 August, 1976 was one which formed in the central 
portion of the northern project area and tracked towards the northwest. The synoptic 
summary is given in Figure 6.19 and the environmental conditions are shown by 
the tephigram in Figure 6.20. The storm was associated with a rapidly moving trough 
which came out of the southeast. Many storms occurred on this day which was 
declared a no-seed day in the North. The hailswath associated with this particular 
storm is shown in Figure 6.21. Six reports of hail were received from this storm 
with grape size hail reported as the maximum size. 


PPI scans of the storm at approximately 15-minute intervals are given 
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SYNOPTIC SUMMARY 


15 AUGUST 1976 


SYNOPSIS: THE MOST SIGNIFICANT SYNOPTIC CHANGES OCCURRING DURING THE 
CURRENT SERIES OF S/W WERE ASSOCIATED WITH THIS THIRD RAPIDLY-MOVING 
TROUGH. RELATIVELY STRONG PVA AND LARGE HEIGHT AND THICKNESS FALLS 
PRECEEDED IT. WITH THE TROUGH MOTION BEING 165 aT 60 KPH. A MAJOR 
SURFACE LOW AND FRONTAL WAVE FORMED OVER MONTANA, THUS ADVECTING 
ADDITIONAL MOISTURE INTO THE ERN PORTION OF THE AHP AREA, 


AHP AREA: SURFACE WINDS LIGHT, NW IN THE MORNING BECOMING E IN THE 
EVENING. WITH SE WINDS ALOFT. THE MAIN STORM DEVELOPMENT WAS WELL AFTER 
MAXIMUM TEMPERATURE TIME, MAINLY OVER THE ERN HALF. WITH TOPS EXCEEDING 
14 KM WHILE MOVING TO THE NW. THIS HIGHLY-ORGANIZED STORM CONTINUED 
OVER NIGHT. PRODUCING THE MOST SERIOUS HAILFALL OF THE SUMMER, 


TROUGH (T) OR 12 HOUR 500 MB CHANGE IN HAIL {NO. HAIL 
RIDGE (R) VORTICITY [HEIGHTS *|THICKNESS*] TEMPS. SIZE REPORTS 
PASSAGE TIME Gice sou) (M) (°C) |CATEGORY| N/ S 
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Fig. 6.19: Map showing the synoptic situation and a brief Synoptic summary 
for the storm on 15 August, 1976. (From Deibert (ed.), 1977). 
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Fig. 16.20: 


STATION: CYRM 
TEPHIGRAM DATE: 16 AUG 76 


TIME: 00002 


Rawinsonde sounding from Rocky Mountain House at 18:00 MDT, 15 
August, 1976. (From Deibert fed.), 1977). 
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Fig. 6.21: Map showing the hailswath resulting from the storm on 
PS. August,, 19/6. 


S4 


4 


vont ie 


y 
it 


hed 
. | : 
' 


ae 


\ pe 
y 
f 
a 


38 


in Figures 6.22 and 6.23. Reflectivities in excess of 50 dBZ were recorded. 

The radar constant used for the day was 74.2 dBm with a peak trans- 
mitted power of 227 kw. 

The 7 time series for the storm is given in Figure 6.24. The dashed 
line represents the trend in the data computed using the Tukey filter described in 
Appendix A. 

As was the case with the 12 August, 1976 storm there was a period of 
time during this day in which the radar was on a full three-minute scan. The 
period was from 18:07 MDT to 19:47 MDT. Thirty two points of a total 86 data 
points for this day were interpolated. This is not a desirable situation. However, 
as mentioned earlier, previous time series possess autocorrelation coefficients 
greater than .90 for 1.5 minute time lags and, therefore, linear interpolation was 
felt to be preferable to determining the power spectrum from data spaced three 
minutes apart. 

The autocorrelation coefficients of the data for lag one (1.5 min) to 
lag ten (15 min), computed using (5.4.2) are given in Table 6.3. 

The unsmoothed power spectrum of the raw data is given in Figure 6.25 
and the smoothed power spectrum, using a three-point running average, is given in 
Figure 6.26. The smoothed power spectrum of the filtered highpass data is given 
in Figure 6.27. The values of the red-noise power are applied to the raw data 
smoothed power spectrum using an autocorrelation value P = .937, and to the high- 
pass data smoothed power spectrum using 9 = .695. There were no statistically 


significant cycles evident at the 95 per cent confidence level assuming a Chi-square 
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PPl echoes for the storm on 15 August, 1976. 
Reflectivity = from 20 dBZ by 10 dBZ. Time 
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Fig. 6.23: PPI echoes for the storm on 15 August, 1976. Elevation ~ 2.0° 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 19:00 - 19:47 MDT. 
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Fig. 6.25: Unsmoothed power spectrum of n for the storm on 15 August, 1976. 
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Fig. 6.26: Three-point Daniell-filtered power spectrum of n for the storm 
on 15 August, 1976. The red-noise spectrum is shown. 
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Fig. 6.27: Three-point Daniell-filtered power spectrum of Tukey-filtered 
n highpass data for the storm on 15 August, 1976. The red-noise 
spectrum and 90% confidence interval are shown. 
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Table 6.3: Lag 1.5 to 15 minute autocorrelation coefficients of n 
for the storm on 15 August, 1976. 


LAG NUMBER TIME LAG (min) AUTOCORRELATION COEFFICIENT 


l [i Sy faye: 
2 3.0 .840 
3 4.5 2701 
4 6.0 2070 
5 Tino 45] 
6 9.0 34] 
7 | 10.5 Paks: 
8 12.0 105 
9 13.5 - 009 
10 15.0 ~.108 


distribution, although a preferred frequency appears to exist in Figures 6.26 and 
6.27 near 0.025 cycles per minute. This corresponds to a cycle period of 40.0 
minutes. 

Red-noise cannot represent the data completely in this case since the 
autocorrelation coefficients were negative at long time lags. By definition 
R(T+ 1) =R(1)° which approaches zero for large lag times 1, but never acquires 
negative values. The red-noise representation is still superior to that of white noise. 
6.5 The 20 August, 1976 Storm 


The storm on 20 August, 1976 formed just inside the northwestern project 


boundary and tracked eastward, dissipating in the northeast quadrant approximately 


2.5 hours later. The storm was seeded. 


ie eae Ge Aaa Gee 
rT v 


Crisis 900 D Metta if 


iN 
i —_ o 
eel) 
; Tals ‘ 
| + 
i ‘sae A 
Uy; 
| 
+ 
: lj 
aN an 
y OOK}, = 
= 
‘ 
/ 
A ma ti 
| ';. r .) i § 6 
J - t arr. ae 
” in ‘i ti nr 5 ¢ ’ Tala Oo 125 
a 4) va ‘ . . | af ' -, ~ a aT 
PHT Par | “i a “prin i‘ fol To evitonon ay pcre zingial lie 
; ’ 7 7 a P ; % is ‘* —, ’ rz j % i >) \¢ ig . gs 
| Vivpop taven ted ,% gamit gms gral vel ces iarlooorgqe si { 
ve i ; - ¥ ; pe ; - : . ee Ai 
+f ; 7 7 % r * : » . ‘-_ ® , oa 
dian Sthiw ie Hore of vermque | iz a mohin inane er G2 bon beg =’ oT eSUF 
7 j i : 7 : i 
i> : ; a: h i > ae mach in, 
: > “3 : , : - meas 1 wh a 7 * " 
Dae mw : sae — : 


1 


,) foals * 4 laa a wor ‘wrth abit pul loeenetey ons yt mabe! no my 


non a 
© yletanixciage wen Atgeiltyo0 ait a othe denen pa f 


ac. 4 
nw) - var aN. . 2 : eee 


Beis 


102 


The synoptic summary is given in Figure 6.28 and the environmental 
conditions are shown Ey the tephigram in Figure 6.29. The resulting hailswath 
is shown in Figure 6.30 which indicates that 8 hail reports were received, all re- 
porting pea-size hail as the maximum size. 

PPI scans of the storm at approximately 15 minute intervals are given 
in Figures 6.31 and 6.32. The storm started as three tiny centers and grew into 
two significant precipitation cores which later joined to form one major storm 
center. Maximum reflectivities were slightly in excess of 40 dBZ. 

The 7) time series for 20 August, 1976 is given in Figure 6.33 which 
also indicates the time and amount of cloud seeding. The dashed line represents 
the trend in the data computed using the Tukey filter described in Appendix A. The 
radar constant used for the day was 72.3 dBm with a peak transmitted power of 
351 kw. 

The radar operated at a 1.5-minute scan period the entire day and 
therefore, no interpolation of data points was necessary. 

The autocorrelation coefficients of the data for lag one (1.5 min) to 
lag ten (15 min) computed using (5.4.2) are given in Table 6.4. The time series 
displayed very high autocorrelation despite the storms rather complex nature. 

The unsmoothed power spectrum of the raw data is given in Figure 6.34 
and the smoothed power spectrum using a three~point running average is given in 
Figure 6.35. The values of the red-noise power are applied to the smoothed power 
spectrum using an autocorrelation value p = .985. 


The highpass data smoothed power spectrum with the red-noise power 


| were srt | 
iG 
ovie om @ vienait deivaovt uletonivonge ta cro of 

OWe wire INTD evi oreapane mn bein mivate on ta baw { 


an i mnie ma bain Ta vats oti aaron rte 


i 1 i 
eu o«@ 


ms w Oe ig ino, toy th bal Bieta meee eaitivitas lien i 


\ 4 et a i. 
whey be oP apie iF ene densi ect xt | { teal ON e i rahe: conti r aa 


soa! bora aal . eilbeos tum tata ail si o sri ata 
en * q ' ' j 1. & { a a ye 
fas c’) ing A a Fa . 1 rs ti) ch. D ! iy Ye MU f ail’ Gere net hare ig a uta allt nt =f 
i rc " . 7 7, _ | . . . ‘ 
loo nwoe athe sooq aiitiw mob 2% gow yop elt to? eee ae 


eS f ft 


=a qemty ovyl eves rtd heey ern ee} vet. | o-38 bsigiajo whet od? 
7 —_ 


; pe : 7 + a 
; | Pp SOMES Keim ig erin stib te coitalempain on anal 


s J i | 7 


‘tdedin BP ete agal ast? otet off Io etnel Hage oeitalamereted oll - 


‘ oa 
‘Soha aloe b4 ote ni nevig ew ie a) pri bonnie (alin 

ewtan es Lyre, prt parse awit gliaweb we colwaaaclve dat ai vy 
bg 2 an Widevigl ue EF ont te mane owen biteaiie aia 


| . owt nag # re gnionn tHe o prize marion owegs 


0 


} 7 yl ing, eadicoms sat 12 aes “19 wed ation att te sl r 


ae! 


a) 
1 


$ ge ‘efter m me 


‘i 


103 


SYNOPTIC SUMMARY 


20 AUGUST 1976 


SYNOPSIS: THE MOST IMPORTANT SYNOPTIC FEATURE HERE WAS A RAPIDLY-MOVING 
TROUGH PASSING THROUGH AHP BEFORE NOON. (NOTE: 12-HR CHANGES SHOWN ARE 
FOR 20/00-12Z.) THE APPARENT WEAK RIDGING BEHIND THE TROUGH WAS QUASI- 
STATIONARY, HENCE NOT A S/W FEATURE. THIS RESULTED FROM THE START OF 
RETROGRESSION FROM E TO W OF THE L/W RIDGE. 


SURFACE PRESSURE FEATURES WERE LITTLE-CHANGED FROM THE ISTH. WITH 
SOUNDINGS FOR BOTH DAYS SHOWING THAT THE INCREASE IN INSTABILITY WAS A 
RESULT OF UPPER-LEVEL CHANGES. 


AHP AREA: W WINDS AT ALL LEVELS. SEVERAL ORGANIZED CELLS FORMED ALONG 
THE FOOTHILLS NEAR NOON, THESE DEVELOPED TOPS TO 9 KM WHILE MOVING 
RAPIDLY E THROUGH THE AREA AND DECAYED AFTER MAXIMUM TEMPERATURE TIME, 


TROUGH (T) OR 12 HOUR 500 MB CHANGE IN HAIL NO. HAIL 


RIDGE (R) VORTICITY |HEIGHTS * {|THICKNESS*} TEMPS. SIZE REPORTS 
PASSAGE TIME |(107°S7!) | (m) (M) (°c) |caTecory| N/ S 


500 MB ANALYSIS _— 
HEIGHTS AND VORTICITY an 
00002 21 AUG 1976 


Fig. 6.28: Map showing the synoptic situation and a brief synoptic summary 
for the storm on 20 August, 1976. (From Deibert (ed.), 1977) 
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Tp STATION: CYRN 
en GRAM ee oa aie 76 
TIME: 00002 
180° 


Fig. 6.29: Rawinsonde sounding from Rocky Mountain House at 18:00 MDT, 
20 August, 1976. (From Deibert (ed.), 1977) 
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Fig. 6.31: PPI echoes for the storm on 20 August, 1976. Elevation ~ 2.0° 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 19:00 - 20:15 MDT. 
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Fig. 6.32: PPI echoes for the storm on 20 August, 1976. Elevation ~ 2.1° 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 20:30 - 21:15 MDT. 
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POWER 
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0.000 90.033 0.066 0.999 0.132 0.165 
CYCLES PER MINUTE 


Fig. 6.34: Unsmoothed power spectrum of n for the storm on 20 August, 1976. 
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Fig. 6.35: Three-point Daniell-filtered power spectrum of n for the storm 
on 20 August, 1976. The red-noise spectrum is shown. 
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Table 6.4: Lag 1.5 to 15 minute autocorrelation coefficients of n 
for the storm on 20 August, 1976. 


LAG NUMBER TIME LAG (min) AUTOCORRELATION COEFFICIENT 


l 35 S773 
2 3.0 Sie 
3 4.5 .898 
4 6.0 5082 
3 TAS) | ./64 
6 8) .690 
7 | 10.5 n022 
8 12.0 048 
2 psi) 473 


10 15.0 Lor? 


applied to it using an autocorrelation value of © = .740 is given in Figure 6.36. 

There were no statistically significant peaks at the 95 per cent con- 
fidence level assuming a Chi-square distribution. However, preferred cycles were 
displayed at 0.031, 0.057, and 0.089 cycles per minute. These frequencies cor- 
respond to cycles with periods of 32.3, 17.5, and 11.2 minutes. 

This storm was seeded with 8 cloud-base flares burned over a period of 
approximately 1.25 hours. 

A feature of this storm which contrasted with the other hailstorms in this 
study was the very slow rate at which it dissipated. No other storm displayed such 


a long drawn out decline in the trend. The storm passed over Pigeon Lake during 
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noise spectrum and 90% confidence interval are shown. 
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Three-point Daniell-filtered power spectrum of Tukey-filtered 
n highpass data for the storm on 20 August, 1976. The red- 
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the period 20:30 MDT to 21:00 MDT but the n time series did not display any 
change during this period which could be credited to a lake effect. 
6.6 The 24 August, 1976 Storm 

A severe storm occurred in the northern area on 24 August, 1976. 
The storm formed in the foothills region near the western boundary of the project 
area and tracked northeastward. The synoptic summary is given in Figure 6.37 
and the environmental conditions are shown by the tephigram in Figure 6.38. This 
storm was not seeded. 

The hailswath for this storm is shown in Figure 6.39. The storm pro- 
duced sporadic patches of hail, indicating that it experienced several cycles in 
intensity. Walnut size hail was reported in two instances, however, the majority 
of the hailswath indicates pea-and grape-size as the maximum sizes. 

PPI scans of the storm at approximately 15-minute intervals are given 
in Figures 6.40 to 6.44. Reflectivities in excess of 50 dBZ were recorded. 

The 1 time series for 24 August, 1976 is given in Figure 6.45. The 
dashed line represents the trend in the data computed using the Tukey filter des- 
cribed in Appendix A. The trend indicates at least three major intensifications in 
7” during the storm lifetime. The radar constant for the day was 74.0 dBm with 
a peak transmitted power of 238 kw. 

The radar operated on a three-minute scan for the period 17:20 MDT 
to 19:08 MDT and consequently 37 data points had to be interpolated. 

The autocorrelation coefficients of the data for lag one (1.5 min) to 


lag ten (15 min), computed using (5.4.2) are given in Table 6.5. 


Oe se i edn 


5 85 enit a no af ere Fane : ) | 
MT 8.0 swith te agin ot om resort oll csi Te 


F 
fi 
ive 
; 5 a 7 
o thy 
ae ‘ zs > 
* i 


: ian. La ae TAL 
meng meinte th <QBid engi Yi wha al miots tie adh Ahowtiod at e 
i | — ; ae 
gl: wee besnetacre UT holt edtmaibel’ sort ioe vetotag a 
visors 2/7 avenond _ CAFS moterti cue§ avi biel wyo¥ ew liot oxia ial 7 : cite 


i 
} 


; ~~ ‘ e . i z 
ANTS MUM KON St so SS eqnip erro ~ 004 ssinsilne ei 
iy > 4 ‘ ® * ; . A - 7 it, - a’ é eo r 
svi 610 stovistal Siuinrets vy atom xoigen-Te rials art ie enna #94 a. : 
: i ae : i, 


; y 2 . 4 
= a : od Cae 
bobroosy eiew SSb 02 to exsoas of poltivitsellon .0h,.aa8 0b, a: ; 


. -; ae : - yt ite oe 
oft be qual fi ctawip #1 SVGT tuguA BR vet astree mitt oF ‘oft 


-eob 10 lit yooluT art ghiew betugmon etoborit ol nari art etnozonqay 


“4 
- 


ni anal fooltjanaitel sojain went tanel to eattal | orate nal ia et it 
Asie infil at , ace - ait 1 inotenoa sober ‘oft | cae fis mote 
my _ oy  awsl eaten | 


TOM ag.Nt boireq ‘ot te mene e stuninesondtc Bre beteweis sobend 


i. kas _hetelagnarat od. of biel aniog piob ct viomceots ba 
ar ey ea ihn. B anit go! otal anh ¥e snaiaitieos sles ad i 


8 fox chia Whew c. +0 eit ° 


15 


SYNOPTIC SUMMARY 


24 AUGUST 1976 


SYNOPSIS: L/W FEATURES, VIZ THE RIDGE OVER MID-CONTINENT AND THE E 
PACIFIC TROUGH, GAINED AMPLITUDE FROM AUGUST 21-23, FORCING S/W FEATURES 
N OF ALBERTA. THE RIDGE STARTED TO COLLAPSE ON THE 24TH, THEN MOVED E. 
ALLOWING THE MAIN 500 MB STREAM TO DROP S OVER CENTRAL ALBERTA WITH AN 
EMBEDDED S/W TROUGH. A SURFACE TROUGH OF LOW PRESSURE FORMED ALONG THE 
FOOTHILLS IN RESPONSE TO THIS S/W, GRADUALLY EVOLVING INTO A LOW CENTER 
BY THE NEXT MORNING. 


AHP AREA: SURFACE WINDS LIGHT S OVER MOST OF THE AREA VEERING TO SW AT 
500 mB. A SINGLE-CELL STORM DEVELOPED SW OF CYRM DURING LATE AFTERNOON 
TRACKING NE 35 KPH, WITH MAXIMUM TOPS 12 KM, THEN DYING LATE IN THE 

EVENING, 


TROUGH (T) OR 12 HOUR 500 MB CHANGE IN HAIL NO. HAIL 


RIDGE (R)  |VORTICITY|HEIGHTS * ee cae SIZE | REPORTS 
PASSAGE Time |(107°s7') | (M) (m) (°c) |caTecorY| N/ Ss 
R-24/19002Z 
T-25/1000Z 


: 


500 MB ANALYSIS 
HEIGHTS AND VORTICITY i 
00002 25 AUG 1976 LAA | 


Fig. 6.37: Map showing the synoptic situation and a brief synoptic 
summary for the storm on 24 August, 1976. (From Deibert 
(eds¥ 1977) 
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STATION: CYRM 
TEPHIGRAM iy. 25 aug 76 


8 Cm 0 TIME: 00002 
%,12 


Rawinsonde sounding from Rocky Mountian House at 18:00 MDT, 
24 August, 1976. (From Deibert (ed.), 1977) 
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Fig. 6.40: PPI echoes for the storm on 24 August, 1976. Elevation ~ 1.9° 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 16:29 - 17:44 MDT. 
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Fig. 6.41: PPI echoes for the storm on 24 August, 1976. Elevation ~ 2.0° 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 17:59 - 18: 44 MDT. 
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Fig. 6.42: PPI echoes for the storm on 24 August, 1976. Elevation ~ 2.1° 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 18:59 - 19:45 MDT. 
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Fig. 6.43: PPI echoes for the storm on 24 August, 1976. Elevation ~2.1° 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 20:00 - 20:45 MDT. 


20:00:22 n = 1.32x10 - cm 20:15:18 Per rinOa Omens 


20: 30:13 n= a AsePalens cm 20:45:09 n= TeRoet OF ween. 


120 


Fig. 6.44: PPI echoes for the storm on 24 August, 1976. Elevation ~ 2.1° 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 21:00 - 21:29 MDT. 
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Table 6.5: Lag 1.5 to 15 minute autocorrelation coefficients of n 
for the storm on 24 August, 1976. 


LAG NUMBER TIME LAG (min) AUTOCORRELATION COEFFICIENT 


1 1 «936 
2 3.0 EOL 
3 4.5 687 
4 6.0 .554 
5 7.9 aay 
6 9.0 396 
, | LOco -365 
8 250 349 
9 13.5 -327 
10 15.0 sol] 


The unsmoothed power spectrum of the raw data is given in Figure 
6.46, and the smoothed power spectrum, using a three-point running average, is 
given in Figure 6.47. The values of the red=noise power are applied to the 
smoothed power spectrum using an autocorrelation value p = .921. The highpass 
data smoothed power spectrum with the red=-noise power applied to it using an 
autocorrelation value of p = .806 is given in Figure 6.48. 

Once again, no statistically significant peaks in the power spectrum 
were found assuming a Chi-square distribution. However, preferred cycles are sug- 
gested by the broad peaks in Figure 6.48 near the frequencies 0.034 and 0.117 
cycles per minute. The corresponding cycle periods are 29.4 and 8.5 minutes 


respectively. These appear to be the major and minor fluctuations in the overall 
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POWER 
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Fig. 6.46: Unsmoothed power spectrum of n for the storm on 24 August, 1976. 
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Fig. 6.47: Three-point Daniell-filtered power spectrum of n for the storm 
on 24 August, 1976. The red-noise spectrum is shown. 
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0.000 0.033 0.066 0.099 0.132 0.165 
CYCLES PER MINUTE 


Fig. 6.48: Three-point Daniell-filtered power spectrum of Tukey-filtered 
n highpass data for the storm of 24 August, 1976. The red- 
noise spectrum and 90% confidence interval are shown. 
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storm intensity which caused the rather sporadic hailfall. 
This concludes the presentation of the August, 1976 storm data. 
Further interpretation and discussion of the data is given in a following chapter 


after all of the storm data have been presented. 


CHAPTER 7 
THESZO2 UIP 7 5GA SES TUDIES 
7.1 Introduction 

Three storms which occurred on 20 July, 1975 were analyzed. All 
three storms formed in the foothills region and tracked eastward across the north- 
ern area. The first storm, designated as storm A, is included in this study for 
comparison purposes since only rain was reported to have fallen from it. The 
second storm (storm B) and the third storm (storm C) are analyzed in order to in- 
vestigate the usefulness of the AHP 1975 digitized data. At times during 1975 
portions of returned echoes are missing due to a computer recording malfunction. 
The result was that individual rays or groups of rays were often lost. An example 
of an echo showing the missing rays is given in Figure 7.1. The gaps did not 
always occur nor did they occur at the same spot in the echoes. Asa result if was im- 
possible to predetermine the error associated with the 1 calculation when missing 
rays were encountered. 

The program used to calculate 7 was modified to take into account 
missing rays. When missing rays were encountered, the bin volumes for the rays 
on adjacent sides, were interpolated across the gap to account for the missing 
volume. The known power values in adjacent rays were assumed to represent the 
increased volume, rather than interpolating the returned power values for the mis- 
sing rays. 

The error in 1 would be minimal if the rays on each side of a gap are 
representative of what was actually there. However, the error is large if the mis- 


sing rays occurred in a region with a high reflectivity gradient. 
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The storms from 1976 displayed 7 time series with lag one autocor- 
relation coefficients greater than 0.90 and no dramatic changes in 7 occurring 
over short periods of time. This characteristic would also be expected of 1975 
storms and, therefore, erratic changes in 7 values would probably indicate the 
misrepresentation of missing data by interpolation. 

7.2 20 July, 1975: StormA 

Storm A is included in this study as a test case. No hail was reported 
to have fallen from it, however, radar reflectivities greater than 40 dBZ were 
recorded during its lifetime. The test was to see if the nN values were sufficiently 
small to characterize a rainstorm as opposed to the 1 values for Hanlerorns already 
discussed. 

PPI scans of the storm at approximately 15-minute intervals are given 
in Figures 7.2 and 7.3. The 7 time series is given in Figure 7.4. The dashed 
line represents the trend in the data computed using the Tukey filter described in 
Appendix A. Times during which missing data was encountered are also indicated. 
The radar constant used for the day was 76.4 dBm with a peak transmitted power of 
137 kw. 

The maximum 7 value obtained was 4.97 x 107% cm=* at 15:39 MDT. 
The cumulative frequency curve for all hail reports given in Figure 4.5 indicates 
that approximately 15 percent of all the hail reports occurred with lower values. 
One may interpret this as suggesting that during the lifetime of storm A the probabi- 
lity of hailfall at the surface was at all times less than approximately 15 percent. 


The autocorrelation coefficients of the data for lag one (1.5 min) to 
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Fig. 7.2: PPI echoes for storm A on 20 July, 1975. Elevation ~ les 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 15:30 - 16:30 MDT. 
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Elevation ~ lees 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 16:45 - 17:15MDT. 
“ym 
\44 
128 is 
=~ 
“ony 16:45:23 
a cy ae) peal eae 
96) 
0 
km 
144 
a, 
128 Sp 
o 
P 17:00:19 
W 
= aa) = 1 
eo n = 2.28x10 Leth a 
96 
rye 


km 
144 
128 
112 & W72t5214 
— - -1 
oy n= 2.86x10 Boe 
96 


Chedard] 


A Ene — = 


ty ar 
an | 
ma ‘hort 
0 a 
a 
a 
Oa 
hc Zfe tt 
| — , te ¥t-, eg 
. a> «=i noe op 
it Oy 
1 , 
ti 
| » 
ii ! 
} ; = 
‘i : A i : Af 7 
yl . Se 7 =. Cala 


os 


131 


*poqzeodipul ove 
pesejzunosue sem ejep bulss}w yosiym Bulanp sow!| ‘*pussz soqzeorpul 


aul, peyseq “G/6, SALNF OZ UO Y W4O IS 4O¥ DWI. SNS4dA LU ey Dis 


(LGW) SWIL 


Glidl Gp.Sl cI.9l GviSi GI.Sl 


"ol —=—— = =e 
(VLVG ONISSIW) 


Oo 
5 "92 
G) 
SI 
44) 
roe | 


V WHOLS (‘GZ6I AIAN Oe 


Att 
“we ; tte amit - ~ = 


e 


a 7 x = + P = 
' Od oi us i oy ¥ 4 
onli beled .<\¢' bobs e723 t sua wha t 5 
- x » . PT ti. me rm | b 7 bina wa Le ok 
beistrpcons esx ota’ orfiazin ao fw G 2 a 
east . bat . = . 


132 


lag ten (15 min) computed using (5.4.2) are given in Table 7.1. 


Table 7.1: Lag 1.5 to 15 minute autocorrelation coefficients of 7 


for Storm A on 20 July, 1975. 


LAG NUMBER TIME LAG (min) AUTOCORRELATION COEFFICIENT 


I lz 404 
2 3.0 -406 
3 4.5 -400 
& 6.0 205 
5) ) Tie 2107 
6 wee .076 
7 . TOPS .021 
8 1250 - .022 
9 V3.0 .001 
10 1520 -.162 


This storm displayed the least autocorrelation at short time lags of 
any of the storms studied. 

The unsmoothed power spectrum of the raw data is given in Figure 7.5 
and the smoothed power spectrum using a three-point running average is given in 
Figure 7.6. The value of the white-noise power is applied to the unsmoothed 
power spectrum and the values of the red-noise power are applied to the smoothed 


power spectrum using an autocorrelation value 9 = .533. The white-noise spectrum 


is included since of all storms this one displayed variations closest to random. The 
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Fig. 7.5: Unsmoothed power spectrum of n for storm A on 20 APE, 28 497 Die 
The white-noise spectrum is shown. 


POWER 


0.000 0.033 0.066 0.098 0.132 0.165 
CYCLES PER MINUTE 


Fig. 7.6: Three-point Daniell-filtered power spectrum of n for storm A 
on 20 July, 1975. The red-noise spectrum and 90% confidence 
interval are shown. 
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red-noise spectrum which accounts for the persistence, reveals no significant 
periods in the caso thedien te: spectrum. 

The low energy storm on 20 August, 1976, which was discussed in 
Section 6.5 and which resulted in 8 reports of pea size hail at the surface, dis- 
played similar values of 7. The maximum 7 value for 20 August, 1976 was 
6.56 x 107** cm=", corresponding to a hailfall probability of slightly greater than 
15 percent. However, the majority of the 7 values were less than those of storm 
A. It appears that the parameter 7 is not sensitive enough to distinguish small 
hail from rain with any accuracy. 

The spectra from the two storms are significantly different. The storm 
on 20 August, 1976 displayed high autocorrelation as opposed to the near random 
nature of 1 in storm A on 20 July, 1975. Although some of the randomness may 
be due to the interpolation of missing rays, it is the lack of a well defined begin- 
ning and end to the 1 time series which reduces the variance and, therefore, results 
in small values for short time lag autocorrelation coefficients and low frequency 
power. 

7:3 20 July; 19752 Storm B 

The hailswath for storm B 20 July, 1975 is shown in Figure 7.7. The 
maximum hailsize reported was grape size and the hailswath showed patches of 
hailfall implying that the storm experienced several cycles in intensity. Also shown 
in Figure 7.7 is the hailswath from storm C on 20 July, 1975 which will be discussed 
in the next section. 


PPI scans for storm B at approximately 15-minute intervals are given 
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in Figures 7.8 through 7.11. The data for storm B suffered extensively from the 
problem of missing rays in the radar echoes. Figure 7.8 points out several PPI 
scans which include missing data. The contour program interpolates directly across 
the gap resulting in the unrealistic straight contour lines. 

The 7 time series for storm B 20 July, 1975 is given in Figure 7.12 
which also indicates the times and amounts of cloud seeding, and the periods of 
missing data. The dashed line represents the trend in the data computed using the 
Tukey filter described in Appendix A. The missing data had the most serious effect 
over the interval from 17:15 MDT to approximately 18:45 MDT as seen by the rather 
erratic fluctuations in 1). 

The autocorrelation coefficients of the data for lag one (1.5 min) to 
lag ten (15 min) computed using (5.4.2) are given in Table 7.2. 

The unsmoothed power spectrum of the raw data is given in Figure 7.13 
and the smoothed power spectrum using a three-point running average is given in 
Figure 7.14. The values of the red-noise power are applied to the smoothed power 
spectrum using an autocorrelation value p = .845. 

The application of the Tukey highpass filter resulted in a data set with 
an autocorrelation coefficient 9 = .023. The smoothed power spectrum for the 
highpass data with white-noise power applied, assuming a Chi-square distribution, 
is given in Figure 7.15. There are 3 peaks in the power spectrum that are significant 
at the 95 percent confidence level with respect to white-noise. The statistically 
significant frequencies are 0.063, 0.128, and 0.148 cycles per minute. These 


correspond to cycle periods of 15.9, 7.8, and 6.8 minutes respectively. 
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Fig. 7.8: PPI! echoes for storm B on 20 July, 1975. Elevation ~1.1° 
Reflectivity = from 20 dBZ to 10 dBZ. Time = 15:30 - 16:45 MDT. 
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Fig. 7.9: “PP «echoes for ‘storm Bron 20 July, 1975. ‘Elevation ale 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 17:00 - 18:14 MDT. 
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Fig. /710:* PPl echoes: forrstornmeb fonn20 Julyy 19755. Elevation = 133% 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 18:29 - 19:01 MDT. 
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hig. /.'l: PPI echoes for storm Bon 20 July, 1975. Elevation ~ le 
Reflectivity = from 20 dBZ by 10 dBZ. Time = 19:14 - 19:59 MDT. 
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Fig. 7.13: Unsmoothed power spectrum of nN for storm) Beon.20 July, 1975. 
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Fig. 7.14: Three-point Daniell-filtered power spectrum of n for storm B 
on 20 July, 1975. The red-noise spectrum is shown. 
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Table 7.2: Lag 1.5 to 15 minute autocorrelation coefficients of 7 


for Storm Bon 20 July, 1975. 


LAG NUMBER TIME LAG (min) AUTOCORRELATION COEFFICIENT 


| has) .828 
2 3.0 743 
3 4.5 ./34 
4 6.0 .668 
5 load .608 
6 | 9.0 .580 
7 1055 .485 
8 P2.0 43] 
9 13.5 .408 
10 oO 34] 


The effect of the missing data can be seen in the 7) series as spikes 
superimposed on the trend in the data. The missing data does not appear to be 
serious enough to mask the trend but considerable high frequency noise seems to 
have been introduced. We may expect approximately 3 spectral lines out of the 
64 lines in Figure 7.15 to reach significance by chance at the 95 percent con- 
fidence level, therefore, the two peaks in the power spectrum at 0.128 and 0.148 
cycles per minute may well have been introduced accidently by the interpolation 


of the missing data. Any cloud-seeding effects would have to be discussed in terms 


of the trend in the data. 
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Fig. 7.15: Three-point Daniell-filtered power spectrum of Tukey-filtered 
n highpass data for storm B on 20 July, 1975. The white-noise 
spectrum and 90% confidence interval are shown. 
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Peae20 uly, 1975: e5torm.c 

The hailswath for storm C 20 July, 1975 is shown in Figure 7.7. As 
was the case with storm B, grape~size hail was the maximum size reported. 

PPI scans for storm C 20 July, 1975 at approximately 15-minute in- 
tervals are given in Figures 7.16 through 7.18. The storm was not a single well- 
defined unit, but consisted of two precipitation cores for a major portion of its 
lifetime. The storm passed directly over Gull Lake at approximately 19:50 MDT 
after which it gradually dissipated. 

The 7 time series for storm C 20 July, 1975, along with the times 
and amount of cloud seeding, are given in Figure 7.19. Once again there were 
echoes which contained missing rays and these periods are also indicated in 
Figure 7.19. The dashed line represents the trend in the data computed using 
the Tukey filter described in Appendix A. The 1 time series appears to go through 
one major cycle which corresponds with the continuous hailswath at the surface. 

The autocorrelation coefficients of the data for lag one (1.5 min) 
to lag ten (15 min) computed using (5.4.2) are given in Table 7.3. 

The autocorrelation coefficients are similar to those found in the 1976 
storms. It appears that interpolation due to missing data in the calculation of 7 
did not result in serious errors. 

The unsmoothed power spectrum of the raw data is given in Figure 7.20 
and the smoothed power spectrum, using a three-point running average, is given 
in Figure 7.21. The values of the red-noise power are applied to the smoothed power 


spectrum using an autocorrelation value © = .948. The values of the red-noise 
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Fig: jsl/:  PP1l @ehoes Tor storm Con 20 July; 1975; Elevation + led 


Reflectivity = from 20 dBZ by 10 dBZ. Time = 19:16 - 20:00 MDT. 
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Fig. 7.18: PP echoes for storm € on 20 July, 1975. 
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Elevation ~ 1.3 
Reflectivity = from 20 dBZ by 10 dBZ. 
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Fig. 7.20: Unsmoothed power spectrum of n for storm C on 20 July, 1975. 
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Fig. 7.21: Three-point Daniell-filtered power spectrum of n for storm C 
on 20 July, 1975. The red-noise spectrum is shown. 
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Table 7.3: Lag 1.5 to 15 minute autocorrelation coefficients of 7 


for Storm C on 20 July, 1975. 


LAG NUMBER TIME LAG (min) AUTOCORRELATION COEFFICIENT 


1 liao e756 
2 3.0 .889 
3 4.5 83] 
4 6.0 .788 
5 Nes, ./33 
6 2.0 .687 
7 10.5 .656 
8 12.0 -627 
9 | 135 cays: 
10 15.0 S007 


power using 2 = .467 are applied to the highpass data smoothed power spectrum 
in Figure 7.22. A statistically significant peak at the 95 percent confidence 
level exists at 0.063 cycles per minute, assuming a Chi-square distribution. This 
corresponds to a cycle period of 15.9 minutes. 
7.5 Summary of 1975 Data 

The missing data in 1975 definitely hampers an analysis of this type 
which deals with low-level echoes, and which is responsive to small regions of high 
reflectivity. From the cases studied, it appears that the interpolation technique 
used in calculating 1 across the gaps in the digitized data seemed to be sufficient- 


ly accurate that the trend in the data was not lost. The linear interpolation did 
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Fig. 7.22: Three-point Daniell-filtered power spectrum of Tukey-filtered 
n highpass data for storm C on 20 July, 1975. The red-noise 
spectrum and 90% confidence interval are shown. 
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introduce high frequency noise into the power spectrum, therefore, little con- 


fidence should be placed in the small scale variations. 
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CHAPTER 8 
SUMMARY 

8.1 Data Summary 

The data for each storm has been presented and in each case the re- 
sulting hailswath was a reasonable facsimile of what might be reproduced from 
the 7 series using the distributions for maximum hailsize versus 1 from Section 4.3. 
A particularly interesting feature was that the discontinuity of the hailswaths in 
most cases corresponded well with the trend variation in the 7 series. This feature 
was not always obvious from the individual PPI echoes. 

The interpretation of the power spectra, calculated for each storm, 
must be done with caution since each storm was slightly different in nature. The 
spectra from the storms which were complex and consisted of two or more precipita- 
tion centers should not be compared directly with those of storms which maintained 
a well-defined identity throughout their lifetime. However, each spectrum does 
reveal the essential character of the "| variation for that particular storm. Although 
peaks in a power spectrum may not reach statistical significance at the desired con- 
fidence level with respect to the red-noise power spectrum, they still suggest the 
possibility of preferred cycles in the data. As an attempt to generalize the power 
spectra data for all storms, a procedure of averaging the highpass data periodograms 
to obtain one spectral estimate of the data was performed. This method is discussed 
by Welch (1967). Averaging was applied only to the 1976 storm cases since it was 
felt that the interpolation of missing data in the 1975 case studies introduced an 
excessive amount of high frequency noise in the power spectrum. The 7 storm spectral 


average for the 1976 case studies is given in Figure 8.1. 
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Fig. 8.1: Three-point Daniell-filtered average power spectrum for the 
seven storms of 1976. The red-noise spectrum and 90% confidence 
interval are shown. 
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In general, a large scale cycle in 7 was indicated with a frequency 
of .031 cycles per et corresponding to a cycle period of approximately 30 
minutes. Secondary preferred frequencies were also found in the data with fre- 
quencies of .065, .091, and .104 cycles per minute. The corresponding cycle 
periods are 15.4, 11.0, and 9.6 minutes, respectively. The presence of these 
peaks suggests several scales of variability. On one hand, 7 seemed to vary with 
the overall evolution of the storms as they went through their life cycles. On the 
other hand, small-scale periodicities in 7 were superimposed on the trends. Since 
1 can be related directly to the precipitation falling from cloud base via the mete- 
orological radar equation, this implies that the type, and/or, amount of precipi- 
tation displayed this pulsating state which was regular enough that it could be dis- 
criminated throughout the lifetimes of all the storms. 

The author believes that these small-scale fluctuations in 1 are related 
to the fine-scale features of storms reported by Barge and Bergwall (1976), who pro~ 
posed a concept of a large circulation cell upon which are superimposed smaller 
perturbations. 


8.2 1, as a Hail Suppression Evaluation Parameter 


Since 7 has been shown to be related to the precipitation at cloud base, 


and appears to be correlated with the precipitation at the surface, an attempt was 
made to use 7] as a measure of cloud-seeding effects. However, the results were 


disappointing. 


Five of the hailstorms analyzed were seeded. The amount of cloud seed- 


ing is indicated within 15 minute intervals on the 1 versus time plots for 17 July, 1976 
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(storm A), 6 August, 1976, 20 August, 1976, 20 July, 1975 (storm B), and 20 
July, 1975 (storm C). In all cases except storm A of 17 July, 1976, the storm 

was studied for an ample period before cloud seeding commenced and it was pos- 
sible to study storm behavior before and after seeding. If cloud seeding altered 

the precipitation mechanism significantly, an effect should be displayed by a 
change in 7. Four hailstorms analyzed were not seeded, therefore, the natural 
variability of 7 could be studied for comparison purposes. Seeding rates and seed- 
ing techniques are not emphasized in this study. 

According to Warner (1976), a decrease in reflectivity might be found 
in intense precipitation cores of individual echoes roughly 25 minutes after seed- 
ing. This estimate assumes that a sufficiently high concentration of embryos develop 
in approximately 8 minutes. Another few minutes are required for the embryos to 
be carried aloft and 10 to 15 minutes are allowed for fallout. Seeding effects are 
sought first in the trend of the data which is related to the overall storm intensity. 

In storm Bon 20 July, 1975 (Figure 7.12), 7 increased (16:15 - 16:45) 
when the seeding commenced and, within 30 minutes, 1 decreased. No seeding 
occurred during the interval 17:15 to 17:30 and, 30 minutes later, 7 is seen to 
increase slightly. Seeding commenced once again, and % decreased and then in- 
creased although not reaching the maximum values achieved earlier. Following 
this, the storm dissipated. The early response of 7) suggests a positive seeding effect 
in this case. 

In storm C of 20 July, 1975 (Figure 7.19), seeding commenced after 


1” had already reached a maximum value. The storm dissipated after 19:30. This 


cannot be attributed solelyto the cloud seeding, since the storm crossed directly 
over Gull Lake at 19:30, and it has been suggested by Warner (1976), that a lake 
may have a dissipating effect on storms. 

Storm A on 17 July, 1976 (Figure 5.9), was a mature storm by the 
time it entered the project area. Cloud seeding and radar detection began at ap- 
proximately the same time. The parameter 7 increased when cloud seeding began 
and, approximately 30 minutes later, 7 values decreased. The storm was heavily 
seeded and 7 decreased steadily for approximately 1.5 hours. At 21:45 71 in- 
creased, then decreased rapidly. The response of 1 approximately 30 minutes 
after the initial seeding once again suggests a positive seeding effect. 

Storm B on \7 July, 1976 (Figure 5.19), was considered to be not seed- 
ed although one cloud-top flare was released into the storm early in its lifetime. 
The 7 time series for storm B did not vary as much as that for storm A and this cor- 
responded well with the continuous hailswath. The dissipation rate was not quite 
as fast for storm Bas for storm A. Since the two storms occurred on the same day, 
developed under similar environmental conditions, and travelled approximately 
across the same region, these two storms constitute the best examples with which 
to look for seeding effects. The possible effect storm A had on the environment of 
storm B would be difficult to assess, and is not considered. The raw data power 
spectra were shown not to be significantly different using the Chi-square goodness- 
of-fit test in Section 5.7. The possibility exists that the greater trend variability 
and sporadic hailswath of storm A were results of the cloud seeding. 


The storm on 6 August, 1976 (Figure 6.5) dissipated approximately 30 
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minutes after the initial cloud seeding, although, the storm was shown to have 
reached its maximum intensity by the time cloud seeding began. Whether or not 
the cloud seeding limited the subsequent maximum intensity cannot be assessed. 
The presence of Buffalo Lake may have contributed to the dissipation of the storm. 

The storm on 20 August, 1976 (Figure 6.33), was also seeded, but, 
as was mentioned in Section 6.5, the storm was not a single entity but consisted 
of two precipitation cores for a major portion of its lifetime. The storm was not 
seeded heavily. It is uncertain whether or not 7 could show a seeding effect in 
such a case since the seeding may affect one part of the storm but be masked by the 
other portion. In any event, the 7 series in Figure 6.33 shows that the storm was 
on a general decline when cloud seeding began. The storm dissipated very slowly. 

The non-seeded storms did not display any consistent traits. In general, 
the non-seeded storms dissipated as rapidly as the seeded storms. The trends in the 
data for the non-seeded storms were not notably different than those for the seeded 
storms. For example, the trend in the data for the storm of 24 August, 1976 which 
was not seeded, was just as variable as the trend in any of the seeded storms. 

The technique of averaging the power spectra for the seed and no-seed 
cases was performed. Figure 8.2 gives the 3-point smoothed average power spectrum 
of the highpass data for the 4 non-seeded storms in 1976. The values of the red-noise 
power were applied to the spectrum using the average autocorrelation value for 
the non-seeded highpass data, 9 = .622. A Chi-square distribution was assumed 
with the number of degrees of freedom equal to the sum of the degrees of freedom 


for the individual samples. Similarly, the 3-point smoothed average power spectrum 
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Fig. 8.2: Three-point Daniell-filtered average power spectrum for the four 
non-seeded storms of 1976. The red-noise spectrum and 90Z 
confidence interval are shown. 
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for the 3 seeded storms in 1976 is given in Figure 8.3. The values of the red-noise 
power, assuming a Chi-square distribution, are applied to the spectrum using the 
average autocorrelation value for the seeded storm highpass data, 9 = .7I1. 

There are no significant peaks in either averaged power spectrum. 
Perhaps the Chi-square distributed red-noise model cannot adequately represent the 
data. Preferred cycles are suggested in the seed-and no-seed averaged spectra. 
Their appearance in many storms adds to their significance. Both the seeded and 
non~seeded average spectra displayed a preferred frequency near .033 cycles per 
minute. The corresponding cycle period is 30 minutes. This peak may have resulted 
from the highpass Tukey filter which has effectively removed all cycles with fre- 
quencies less than .033 cycles per minute. The average spectrum for the seeded 
storms displays a well-defined peak at .065 cycles per minute (cycle period: 15.4 
minutes), and a minor peak at .091 cycles per minute (cycle period: 11.0 minutes). 
The average power spectrum for the non-seeded storms indicates a broad, poorly 
defined peak near .055 cycles per minute (cycle period: 18.2 minutes) and a well 
defined peak at .105 cycles per minute (cycle period: 9.5 minutes). 

A distinction between the two averaged spectra is exhibited by the 
presence of substantial power near .105 cycles per minute in the non-seeded cases, 
and the absence thereof in the seeded cases. A hypothetical physical reason for this 
feature, if it was real, may be that cloud seeding was able to upset the regular 
production of high reflectivity pulses in the early stages of the developing clouds 
as they entered the main storm. However, in general, the comment by Panofsky 


and Brier (1963) that, "Cloud seeders claim they can change the spectrum significantly", 


was not demonstrated by the parameter 1 in this study. 
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3: Three-point Daniell-filtered average power spectrum for the 
three seeded storms of 1976. The red-noise spectrum and 90% 
confidence interval are shown. 
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| CHAPTER 9 

. CONCLUSIONS AND RECOMMENDATIONS 
9.1 Concluding Remarks 

This study has attempted to develop and test a radar parameter which 
is able to depict storm intensity and precipitation at the surface better than radar 
reflectivity alone, with hopes that it may be used to detect and measure cloud- 
seeding effects. 

The parameter developed was volume-averaged radar reflectivity 7. 
lts application to low-level PPI scans of 9 hailstorms yielded the relationship 
between 7 and maximum size hail at the surface that was described in Section 4.3. 
The cumulative distribution curves for Naccording to maximum hail size reported 
at the surface in Figure 4.6, allow one to associate a probability of occurrence 
of hailfall of a particular maximum size with nvalues. This removes the depen- 
dence on hail reports and hailpads as the only source with which to define a 
hailswath. The n-hailsize relation should be tested in real time as an instanta- 
neous indicator of what is happening at the surface. 

Since 1 and radar reflectivity factor are related, some of the problems 
associated with their use in hailstorms are comparable. It was possible to distinguish 
between rain and large hail using 7; however, 7 values were sometimes similar 
for rain and rain mixed with small hail. The Nvalues were measured in this study 
using a 10.4 cm wavelength radar and were found to increase monotonically with 
size for hail less than approximately 3 cm in diameter. This monotonic increase 
was no longer evident for hail larger than 3 cm in diameter. The backscattering 


from hydrometeors which are no longer small with respect to the wavelength of the 
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radar exhibits oscillations with increasing particle diameter, therefore, the 
measured 1) values for hail were consistent with theoretical calculations. 

The question arises as to whether a decrease, due to cloud-seeding, 
in the volume of all contours corresponding to "mostly hail", may be masked by 
the concomitant increase in the volumes of all contours corresponding to "mostly 
rain". If the difference in radar reflectivity between "mostly rain", and "mostly 
hail" regions is at least 10 dBZ, the overall volume of the "mostly rain" contours 
would have to increase by at least a factor of two for no net effect on nN. The 
probability of this occurring is not high. Future work for further refinement of 
this technique should involve the calculation of 7 for only those volumes in the 
storm echo which exceed some assumed threshold radar reflectivity value for 
"mostly hail"; for example, 40 dBZ. 

The 1975 radar data was found to be of sufficient quality for a study 
of this type if only ae trend in the data is required. The small-scale or high 
frequency fluctuations in 7) were greatly affected by the need to interpolate 
across missing portions of the data. 

The volume averaging technique allowed storms of different sizes and 
at different ranges from the radar to be compared on even terms. peariaee the 
greatest consequence of this study is the documentation of the variations in 7 
with time. A spectral analysis of the 7 time series for the storms revealed pos- 
sibilities of preferred cycles. Fluctuations in 7, a parameter related to the 
instantaneous precipitation from the storm, were shown to imply a storm concept 


consisting of a large circulation cell with periodicity of approximately 30 minutes, 


upon which is superimposed small-scale perturbations. The small-scale perturbations 
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have periods of approximately 10 and 15 minutes. Attempts to find statistically 
significant differences between seeded and non-seeded cases were hindered by the 
small sample size. Future work should involve the analysis of many more seeded 
and non-seeded storms. 

This study was limited to well-defined storms in order to formulate an 
accurate 7~versus-hailfall relation. However, the technique can be applied to 
all situations. An 7 value can be computed for particular portions of massive storm 
systems by merely defining the "window" to include only those areas of interest. 
Changes in 7 with time may be more meaningful than the absolute value, allowing 
one to monitor storm intensification or dissipation with time. 

In most cases 7 was shown to be more sensitive to overall storm inten- 
sity than radar Ha alone for revealing trends in the data not obvious from 
the PPI scans. The analysis of the changes in 7 with respect to cloud-seeding 


events has been discussed. Unequivocal cloud-seeding effects were not detected. 
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APPENDIX A 
THE TUKEY LOW-PASS FILTER 

If an original time series contains frequencies or periods which are 
of no interest, these waves may be reduced or eliminated by statistical filtering. 
Therefore, if the smooth underlying trend of a time series is desired, the data 
may be filteredto eliminate theshort period variations. A filter of this type is 
called a low-pass filter since the low-frequency or long-period waves are barely 
affected. 

It is possible to design statistical filters such that only frequencies 
greater than a specified frequency are eliminated and the amplitudes of all waves 
in the time series affected by the filtering are known. This information is repre- 
sented by the frequency response function of the filter. 

The smoothing function used in this study is known as the Tukey low- 
pass filter. The filter consists of a series of weights w(t) which are cumulatively 
multiplied by consecutive values of the time series to obtain the filtered variable. 


The weight function is given by (A.1.1) and is illustrated in Figure A.1. 


Tt 
w(t) = he (1 + cos =) | t| <M (A.1.1) 


it} >M 


for lag time t up to a maximum time lag M. The frequency response of the weight 


function is obtained by taking the Fourier transform (A.1.2). 


i -i2 ft 
W (f) || Wit)e dt (AST?) 


am CO 


The frequency response W(f) is given in (A.1.3). 
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Fig. A.J: The Tukey filter weight function. 


Ww (f) 


Fig. A.2: The Tukey filter response function. 
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The frequency response is plotted in Figure A.2. 

The trends in this study were computed using M = 15 minutes. The 
filtered variable was determined using 10 data points on each side therefore 10 
data points were lost at the beginning and end of each low-pass filtered series. 
Thus all frequencies greater than 1/15 = .0667 cycles per minute were eliminated. 
_ The trends displayed in the original time series consist of waves with frequencies 


less than .0667 cycles per minute. 
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APPENDIX B 
THE POWER SPECTRUM OF RED-NOISE AND WHITE-NOISE 
A power spectrum which displays an even distribution of relative 
variance across all frequencies is known as a white-noise spectrum. 
According to Kanasewich (1975), a time series is called white-noise 


if it satisfies the following criteria: 


1. The time average of the series is zero. 
N 
E(x,) = lim Ss > 4 = 0) 
eloha ia Wuaeranuh 


Where E(x, ) is the mean value. 


2. The autocorrelation with zero lag is finite. 
A_(0) = lim 
x t —>0 


@ is the variance of the time series and also the power 


The quantity o 

of the white-noise. 
3. The autocorrelation of the time series is zero for all lags different from 

zero. 

A random time series is an example of white-noise. 

The power spectrum of red-noise has been described by Gilman et al. 
(1963). Specifically it is a first-order linear Markov process. This is a type of 
random process in which the future development depends upon the most recent state. 


The variate x(t) is given by: 


x(t) =a x(t-I)te, 
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where 1 is the autocorrelation at lag one and €4 is some random input at time 
t. The red-noise spectrum is given by: 


L='c8 


oat) epee ee eee 
1 +9°- 2 acos (2Tf) 


The striking characteristic of the red-noise power spectrum is the general suppres- 
sion of relative variance at high frequencies and consequent inflation at low fre- 
quencies. 

It was found that the spectra of 1 for the storms studied fit the red- 


noise spectral model better than the white-noise spectral model. 
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APPENDIX C 


COMPUTER PROGRAM 


The following is a listing of the computer program used 


in. the caleulatton of n. 


INTEGER TIME(2,3),WINDOW( 13), WINDOX(13),AZ(3),EL(3),HDR(41) 
*,DIFC147),CTIME(4),TEMP,A,B,C,EOW, EOF 

REAL TABLE(147),TABLE2( 147) 

LOGICAL*1 RAY(2,147,3),NAME(10),DATE(9),YR(2) 

EQUIVALENCE (YR,1) 


FORMAT(C'’SRADAR CONSTANT? ' } 
FORMAT(F10.0) 
FORMAT(1X,4A2,F8.1,F12.8) 
FORMAT('STAPE DRIVE? ') 
FORMAT(Q,10A1) 
FORMATC'SFILE NUMBER? ‘9 
FORMAT( 110) 
FORMAT(' DATE DATA LOGGED? DD-MMM-YY'/ 
$'S ") 

FORMAT(9A1 ) 

0 FORMAT(///1X,9A1,4X,F10.2) 


m0 Ona TIHEaANe 


Cea a GET RADAR CONSTANT 
WRITE C544) 
READ(5,2)RADCON 
PTH10.**(i227—-RADCON/10.)/(480.*1.75E-6} 
WRITE(5,8) . 
READ(5,9)DATE 

C==-<= COMPUTE TABLES 
DO 50 1=1,147 
R=IRCON(1)/200.0 
TABLE2(1)=R*R 
TABLEC I )=TABLE2(1I)*0.01832 
R=RADCON+20.0*ALOGIO(R) 
DIFC To =Re2 70-200.04+0.5 


D WRIGEGSS Li anABLEC I), TABLECE),R DIF CI ) 
11 FORMAT GING! 7 SE20.655 17.) 
50 CONTINUE 
C<-—=~— ORENS DISK FILE 
CALL ASSIGN(3,'TERRY.DAT' } 
Coe =— GET NAME OF TAPE DRIVE 


WRITE(5,4) 
READ(5,5)LEN,(NAME(I),1=1,LEN) 
CALL SETIDV(NAME,LEN) 

Ca GET FILE NUMBER 
WRITE (5767) 
READ(5,7)I 
CALL FILE(I,HDR,IEND) 

Orde PRINT PARAMETERS 
WRITE(5,10)DATE,RADCON 


C----- INITIALIZE BUFFER POINTERS 
A=1 
B=2 
C=3 

C----- GET FIRST RAY 


CALL XRAY(TIME(1,C),AZ(C),EL(C},RAY(1,1,C0),EOW, EOF ,WINDOW) 
GOTO500 
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C----- ROTATE BUFFER POINTERS AND GET NEXT RAY 
100 TEMP=A 

A=B 

B=C_ 

C=TEMP 


CALL XRAYCTIME(1,C),AZ(C), EL(C),RAY(1,1,€),EOW, EOF ,WINDOW) 


PECEOR SEQ s1 > OR me OWee a) GOTO/00 


i a aati STORE WINDOW INFORMATION 
DO. 120 ies 
110 WINDOXC I )=WINDOWCTI ) 
Cae aee COMPUTE WIDTH OF RAY 
ANGLE=ICLEAR(AZ(B)-AZ(A) )/2+ICLEARCAZ(C )-AZ{B))/2 
Cs<<s= PREPARE TO PROCESS CURRENT RAY 
200 IBINI=WINDOX(5) 
IBIN2=WINDOX( 6) 
D WRITE(5,12)ANGLE 
be FORMATC1X,10F10.4) 
D WRITE(S,13)IBIN1,IBINZ2,WINDOX 
13 FORMAT(1X,1616) 
IF(WINDOX(7).EQ.0}G0TO0220 
Ces-== CONVERT REFLECTIVITY BACK INTO POWER RECEIVED (DBM) 
DOU 220) basset 4i7. 
210 IFCRAY(2,1,8).GT.-100)RAV(2,1,B)=RAY(2,1,B)-DIF{(1) 
Ce-<== INCREMENT SUMS FOR EACH BIN INSIDE THE WINDOW 
220 DO 300 IBIN=IBIN1,IBIN2 


N=RAY(2,1BIN,B) 
IF(N.LT.-90)GOTO300 
POWER=10.0**((N-100)/20.0) 
VOLUME=TABLECIBIN)*ANGLE 
TOTVOL=TOTVOL+VOLUME 
SUM=SUM+VOLUME*TABLE2(IBIN)*POWER 


D WRITE(5,14)1IBIN,N,POWER, VOLUME ,TOTVOL,SUM 
14 FORMAT(1X,217,4E20.6) 
300 CONTINUE 


IF( EOF .EQ.1)GOTO400 
IF(EOW.EQ.0)GOTO100 

C----- COMPUTE AND OUTPUT PARAMETER 

400 IF({ TOTVOL.EQ.0.0)GOTO450 
PARAM=(7.068E+6/PT)*SUM/TOTVOL 


CALL TIMBINCCTIME,TIME(1,B),TIME(2,B)) 


ELEV#ELGBD A107; 
WRITE(3,3)CTIME,ELEV,PARAM 


450 TFCECH EQS I) SOP 

Cse5=— INITIALIZE SUMS 

500 SUM=0.0 
TOTVOL=0.0 

Cse=—< STORE WINDOW INFORMATION 

600 DO-620 7=),13 

620 WINDOX( I )=WINDOWC( I) 

Coa => ROTATE BUFFER POINTERS AND GET NEXT RAY 
TEMP=A 
A=B 
B=C 
C=TEMP 
CALL XRAY(TIME(1,C),AZ(C),EL(C),RAY(1,1,C),EOW, EOF ,WINDOW) 

C=-=== COMPUTE WIDTH OF RAY (NO PREVIOUS RAY} 
ANGLE=23+ICLEAR(AZ(C)-AZ(B)}/2 
GOTO200 

C-=-== COMPUTE WIDTH OF RAY (NO NEXT RAY) 

700 ANGLE=ICLEAR(AZ(B)}-AZ( A) )/24+23 
GOTO200 


END 


176 


PwORH ENAMORED tert vag, 


676 
store 
al Men} 
| - errr mS 
SAET ESA “(90% SAMAMAIDTOSAT LAYS AMC BRA BA 2 “4 a 
' Pi TAR THERES ae 7 o enti. 
a vi ra320an “18 7OS 
rt se a 1nOGHE 2 i} q | wey 


4 a ; ‘at ae (MACE, a3 
| ¥ (b, OEE, kivTAna 
uy OOK IW eSNERT, i simNCel, pf ek 


r (ateat, ¥ vane 
: : “ast yr tint ona. erage ns 
(MOG) DIVTIDTA Nev +4 £ 4948 Reh iF 7m i 
, bah eS a 
cr yaio~ee,€. 594A L SSY AERO R-, *#. : fi. SWASIM 
UoGvid sey wren, HPO WOKS BED 2) U2 THERSAOKNT 
SnIprpanterented oo€ Oe - 


| (avi at, SPvaAse ‘ 
an nore a3 Th, wea n 
(@, OFA ferko > oe, O1eREYO9) 
pOUAM CMTE See PEAT aM IOV Ms 
vl V+ 1007 OT « JOVFOT 
. 23 ey ATO1 907 RAT ReOV Meee 
“MER AOVIOT , ANNIGY, BIOTA Rit). @)7PTAW 
2) rHyOEst, MOS MLE thon > 
fe 7 ; r iu ff oe 
: ; soMOTORC | .03. Kia 
ng leToac6 oy as 
AY MaAhAS THATUO GAA aT BANOS: -- 
HLESTOOEG 0.03, JOVTOV INE - 


= PTOT AUS *tTGXa «3 ate Ta “in 
oH oo {ia aT fy BE METS iY 1 J 7 
Pees 
: SAA 5 AB3,3 re ¢..3Ts 
7 ny oe) #3. roartt: 
i “ a ene. 72) oa 
| ah ; 
- i iP ey j y sli ey 
? | >. Sig Praniow: WORM 
; i : ‘Et, te? Osa 
yf (7 1WO4 Nea ata : 
VAA TREN THe on PASTAS Lie ros 
i - re a * : > sf . 


7 | ae ed J 
(WOCHIW, 308 sna YYAS ene ae ssa 


ENA. WOTV SRS @ ; 
: Sas ea 


a 


ne 
= 


i) ae a - 


gh) > Aes at 
Ds ae ened 
Tete ee We e 


BSS oiea 


Anema tae 


Means Airey 
PY tbe Rh 


We eter a 
Nera 
sy 


rp Sea 
rest pans a 
SMT Ad 
: ey. 
4 hs 
Thera some 
Se ise 


apaie 
44 3: unaabnis 
34% prlinua ty aami 
AAA) A Se om Sm, ye 


